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We extend previous results on the extraction of the Darwin potential to 
all orders in c -2 from the radiation gauge Lienard-Wiechert solution for the 
system of N positive-energy scalar particles plus the electromagnetic field 
in the Wigner-covariant rest-frame instant form of dynamics to the case of 
N positive-energy spinning particles. This is done in the semi-classical ap- 
proximation of using Grassmann-valued electric charges for regularizing the 
Coulomb self-energies and extracting the unique semi- classical action-at-a- 
distance interaction hidden in any Green function used for the Lienard- 
Wiechert solution. By describing semi-classically also the spin of the particles 
with Grassmann- variables, by means of a semi-classical Foldy-Wouthuysen 
transformation applied the the Dirac-like constraints of the manifestly Lorentz 
covariant spinning particles, we determine the coupling of positive-energy 
spinning particles to the electric field in the semi-classical approximation. 
Then we follow the same procedure developed for scalar particles and, in the 
sector where there is no in-radiation, we determine the effective semi-classical 
interparticle potential. Besides the relativistic Darwin term there are spin- 
orbit and spin-spin terms in the potential. Quantization of the lowest order 
(in c -2 ) part of the closed form of the effective Hamiltonian in the case N=2 
reproduces exactly the standard result of the reduction of the Bethe-Salpeter 
equation for the bound states of two spin 1/2 constituents of arbitrary mass 
(hydrogen atom, positronium, muonium). 
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I. INTRODUCTION. 



At the classical level, the accepted mathematical description of the four basic interac- 
tions starts with an action principle. Manifest Lorentz invariance and local gauge invariance 
force these interactions to make use of singular Lagrangians. This implies that their Hamil- 
tonian description must rely on the Dirac-Bergmann theory of constraints . Constraint 
dynamics provides a natural formalism for any program of a unified description of these 
interactions in terms of Dirac-Bergmann observables. Those observables are gauge invari- 
ant and deterministic variables which describe a canonical basis of measurable quantities. 
When one begins with a singular Lagrangian, the canonical momenta pi = dL/dqi are 
not independent. The relations among them, 4> a {q,p) ~ 0, are called primary constraints 
(the weak inequality ~ means that the equality sign cannot be used inside Poisson brack- 
ets). The canonical Hamiltonian H c (q,p) has to be replaced by the Dirac Hamiltonian 
Hp — H c + J2a ^a{t)4> a which, by way of the arbitrary Dirac multipliers, then accounts for 
the restriction to the submanifold defined by the constraints. The time constancy of the 
primary constraints, d t (p a = {0 a ,if£>} « 0, either produces secondary Hamiltonian con- 
straints or determines some of the Lagrange multipliers. This procedure is repeated for 
the secondary constraints and again if necessary (this is the Dirac-Bergmann algorithm) || . 
When the procedure is finished, there is a final set of constraints, Xa ~ 0, defining the final 
submanifold on which the dynamics is consistently restricted, and a final Dirac Hamiltonian 
with a reduced set of arbitrary Dirac multipliers describing the remaining arbitrariness of 
the time evolution. Dirac divides the constraints into two subgroups: i) the first class 
constraints Xa ~ having weakly zero Poisson brackets with all constraints and being the 
generators of the Hamiltonian gauge transformations of the theory and ii) the second class 
constraints Xa ~ Q w hh det({x^, XhlY) 7^ 0' corresponding to pairs of inessential vari- 
ables that can be eliminated. The vehicle that Dirac invented to eliminate, in a systematic 
and symmetrical way, from the dynamical equations those variables constrained by the sec- 
ond class constraints are known as Dirac brackets. Second class constraints may be primary 
or secondary constraints but they can also be introduced by hand. This can be done in 
two different ways. Firstly one can explicitly break the gauge freedom by the introduction 
of gauge fixing conditions on the gauge variables, so that the original first class constraints 
and the gauge fixing ones become a system of second class constraints. Secondly one can 
introduce pairs of solutions to the equations of motions ( for example, field equations whose 
solutions express fields in terms of canonical variables for particles serving as sources for the 
fields) . In general this procedure may break the original manifest Lorentz covariance of the 
theory ( but maintaining Wigner covariance, when we work in the rest-frame instant form 
of dynamics M ) . 

If the end aim is a unified description of the four interactions, then in the application 
of the Dirac-Bergmann algorithm to special relativistic theories it is advantageous to refor- 
mulate the theory in such a way as to allow a natural transition to the coupling of gravity. 
With this in mind, recent studies (see and references cited therein) have taken advan- 
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tage of work by Dirac |JJ on general relativity arriving to a description of special relativistic 
systems on arbitrary spacelike hypersurfaces (parametrized Minkowski theories). Then, for 
every configuration of the system with timelike total 4-momentum the description can be 
restricted to the so-called Wigner hyperplanes orthogonal to the total 4-momentum. This 
is the intrinsic rest frame of the configuration. On it a new instant form of dynamics (the 
Wigner- covariant rest-frame instant form) can be introduced. 

In a recent paper []5| we presented the technical completion of one segment of this pro- 
gram ||, that of N scalar charged particles plus the electromagnetic field when the charges 
of the particles are described by bilinears in Grassmann variables. There we analyzed 
how to extract the action-at-a-distance interparticle potential hidden in the semi-classical 
Lienard-Wiechert solution of the electromagnetic field equations, a subset of the solutions of 
the equations of motion for the isolated system formed by N positive-energy scalar charged 
particles plus the electromagnetic field. The problem is formulated in the Wigner- covariant 
rest-frame instant form of dynamics which is defined on the Wigner hyperplanes or- 
thogonal to the total time-like four-momentum of the isolated system and which requires 
the choice of the sign of the energy of the particles (we considered only positive energies) 

This was possible due to the semiclassical approximation of using Grassmann-valued 
electric charges (Qf = 0, QiQj ^ for i ^ j) as an alternative to the extended electron 
models used for the regularization of the Coulomb self energies. How this happens was 
shown in Ref. |4j], where the Coulomb potential was extracted from the electromagnetic 
potential by making the canonical reduction of the electromagnetic gauge freedom via the 
Shanmugadhasan canonical transformation. This is equivalent to the use of a Wigner- 
covariant radiation (or Coulomb) gauge in the rest-frame instant form. 

Ref. presented the retarded Lienard-Wiechert solution for the transverse electro- 
magnetic field in the rest-frame instant form radiation gauge: in this gauge, due to the 
transversality, the retarded Lienard-Wiechert potential associated with each charged par- 
ticle depends on the whole past history of the other particles. At the semi-classical level 
a single accelerated charged positive-energy particle with Grassmann-valued electric charge 
does not radiate even if it has a non-trivial Lienard-Wiechert potential, avoiding therefore 
the acausal features of the Abraham-Lorentz-Dirac equations, and has no mass renormal- 
ization. However, a system of N charged particles produces, by virtue of the interference 
terms from the various retarded Lienard-Wiechert potentials of the particles, a radiation 
which reproduces the standard Larmor expression for radiation in the wave zone, when the 
particles are considered as external sources of the electromagnetic field and their equations 
of motion are not used. 

If instead the particles are considered dynamical, the use of their equations of motion and 
of the semi-classical approximation leads to a drastic simplification of the Lienard-Wiechert 
potentials and fields. Indeed, if we make an equal time expansion of the delay by expressing 
these potentials and fields in terms of particle coordinates, velocities and accelerations of 
every order, it turns out that all the accelerations decouple at the semi-classical level due to 
the particle equations of motion. Therefore, at the semi-classical level the retarded, advanced 
and symmetric Lienard-Wiechert potentials and electric and magnetic fields coincide and 
depend only on the positions and velocities of the particles, so that we can find their phase 
space expression in terms of particle positions and momenta. 

In this way the semi-classical Lienard-Wiechert potential and fields can be reinterpreted 
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as scalar and vector interparticle instantaneous action- at- a- distance potentials. It is then 
possible to identify a semiclassical reduced phase space containing only particles by elimi- 
nating the electromagnetic field by adding by hand second class constraints which force the 
transverse potential and electric field canonical variables to coincide with the semi-classical 
Lienard-Wiechert ones in the absence of incoming radiation: A± (r, cr) — A±Lw(r,a) ~ 
0, 7t_i_(t, a) — tt±lw( t ,&) ~ 0. Let us remark that this could be done also in pres- 
ence of an arbitrary incoming radiation A x(rad) (r, a) , n ±irad) (r,a) = —j^A^^r, a) 
by modifying the constraints to the form A^{r,a) — A± lw (t, a) — v4^( rad )(r, a) « 0, 
7Tl(t, a) - j?xLw{i~, o) ~ n ±{rad) (r, a) w 0. 

The reduced phase space is obtained by means of the introduction of the Dirac brackets 
associated with these second class constraints. Since the old particle positions and momenta 
are no longer canonical in this reduced phase space, we had to find the new (Darboux) 
basis of particle canonical variables. The generators of the internal Poincare' group inside 
the Wigner hyperplanes in the rest-frame instant form of dynamics can be reexpressed in 
terms of these new variables: the 3-momentum Vn n %\ and the angular momentum <Ju n t) 
become equal to those for N free scalar particles (as expected in an instant form). The 
interaction dependent boosts Ku n t) are proportional to the internal canonical center of mass 
q + inside the Wigner hyperplane: q + ~ are the gauge-fixings to be added to the rest-frame 
conditions V(i n t) ~ , if one wishes to re-express the dynamics only in terms of particle 
internal relative variables. Therefore we have a unfaithful representation of the internal 
Poincare group. Also the energy- momentum tensor has been evaluated in the new canonical 
variables and there is a suggestion on how to find the M0ller center of energy of a cluster of 
n particles contained in the N particle isolated system. 

The Hamiltonian in the rest frame instant form, generating the evolution in the rest-frame 
time of the decoupled external canonical center of mass, is the internal energy generator M = 
VL n1 \ (the invariant mass of the isolated N particle system). The semi-classical Lienard- 
Wiechert solution implies the existence of interparticle action-at-a-distance potentials of two 
types: vector potentials, minimally coupled to the Wigner spin 1 particle three-momentum, 
under the square root associated with the kinetic energies; ii) a scalar potential (including the 
Coulomb potential) outside the square roots. In the semi- classical approximation all these 
potentials can be replaced by a unique scalar potential, which is the sum of the Coulomb 
potential and of a generalized Darwin one for arbitrary N. It is the (semi-classical) non-static 
complete potential corresponding to the one photon exchange tree Feynman diagrams of 
scalar electrodynamics and is a completely new result. The expression we find contains no 
iV-body forces, being simply a sum of two particle interactions. This is a consequence of 
our use of Grassmann charges and of the equal time description on the Wigner hyperplanes 
of the rest-frame instant form. 

In the N = 2 case we obtain a closed form of the solution by evaluating it in the rest 
frame after the gauge fixing q + « 0: the lowest order in 1/c 2 contribution of the generalized 
Darwin potential agrees with the expression of the standard Darwin potential. We then 



2 It is an arbitrary solution of the homogeneous wave equation and must not be interpreted as a 
pair of canonical variables. 
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show that in a semi-classical sense a special solution of the Hamilton equations is the Schild 
solution [H in which the two particles move in concentric circular orbits. 

In this paper we will extend these results to positive-energy charged spinning particles, 
with both the electric charges and the spins described by Grassmann variables (pseudo- or 
semi-classical approximation), whose rest-frame instant form description was given in Ref. 

The resulting generalized Darwin potential will be the (semiclassical) static and non- 
static complete potential corresponding to the one-photon exchange tree Feynman diagrams 
of spinor electrodynamics, after the restriction to the subsector of positive-energy fermions. 

The fact that the positive-energy spinning particles are a pseudo-classical description of 
positive-energy fermions immediately forces us to face a problem absent with scalar particles. 
In the manifestly Lorentz covariant approach, relativistic scalar particles are described by 
means of the mass-shell first class constraints 

Xi = v\ ~ m] « 0. (1.1) 

The positive- and negative-energy sheets of the mass hyperboloid correspond to the particle 
energies p\ ~ ±\Jmf + p 4 ?. When an external electromagnetic field is present we have 

Xi = \Pi-QiMxi)?-™%aO, (1.2) 

so that p® ~ QiA°(xi) ± \Jmf + [ p { — QiA(xi)] 2 . We see that at the classical level the 
positive- and negative-energy sheets of the mass hyperboloid of scalar particles never inter- 
sect even in presence of external electromagnetic fields []. Therefore, there is no problem in 
the rest-frame description of scalar particles with a fixed sign of the energy, since these par- 
ticles will be always interpreted as the classical basis for Tomonaga-Schwinger asymptotic 
states. 

On the other hand, the manifestly Lorentz covariant description of spinning particles 
requires the introduction of extra Dirac-like first class pseudo-classical constraints [ljj 

XDi = Pipti ~ m& 5 « 0, (1.3) 

whose quantization generates the one-particle Dirac equation Q (for each particle we have 
&fi 75, £f 757 M , XDi ~ i-> l^Pi^Y - mi)i/j(pi) = 0). 



3 As shown in Ref. [|T^], Appendix D, scalar particles give the semiclassical description only of those 
solutions of the Klein-Gordon equation coupled to an external electromagnetic field for which 
the Feshbach-Villar approach to the Klein-Gordon field allows a separation of positive-energies 
from negative-energies in an eikonal approximation. Otherwise the non-diagonalizability of the 
Feshbach-Villar Hamiltonian implies that generic solutions cannot be described only in terms of 
scalar particles when external fields are present, but require more elaborate quasi-particle concepts. 

4 The Grassmann variables ^5 describe the spin structure and, after quantization, generate the 
Clifford algebra of Dirac matrices. We assume that the Grassmann variables describing the spin 
of different particles commute: £f£J = for i ^ j, + £f£f = and so on. 
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The mass-shell constraints x% — Pi ~ m ? ~ must be consistent with the pseudo-classical 
version {xm,XDj} — i$ijXi • m this sense the Dirac equation is the square root of the Klein- 
Gordon one. 

When there is an external electromagnetic field, the Dirac-like constraints become 

Xm = \pi? - QiA^Xi)}^ - mite « 0, (1.4) 

and the resulting mass-shell constraints 

Xi = -i{Xix, XDi} = [Pi - QiA(xi)} 2 -m\ + iQiF^xJgg « 0, (1.5) 

contain a non-minimal term connected with the spin-magnetic field coupling in the rest 
frame. 

Since we know that we cannot separate positive from negative energies in the Dirac equa- 
tion coupled to an external electromagnetic field with an exact Foldy-Wouthuysen transfor- 
mation 0, we expect to have the same problems present also to the pseudo-classical level, 
where the pseudo-classical Foldy-Wouthuysen transformation is known [jDJ in absence of 
external electric fields []. Now in the present case this pseudo-classical transformation has 
to be reformulated by using Grassmann-valued electric charges: this further semiclassical 
approximation will allow to get an exact semi-classical separation of the positive and nega- 
tive energies. Therefore, due to spin, differently from the case of scalar particles, already at 
the semi-classical level we must make a pseudo-classical Foldy-Wouthuysen transformation 
which puts the Dirac-like constraint at order Qi in a form producing only even operators 
after quantization. In this way we shall identify the kind of semiclassical couplings of the 
spinning particle to an external electric field which are compatible with the separation of 
positive energies from negative ones. These couplings will appear in the semiclassical ap- 
proximation of the solutions of the one-particle Dirac equation in all those cases in which 
an eikonal approximation plus a standard Foldy-Wouthuysen transformation will allow us 
to identify positive (or negative) energy spinning particle as a realistic approximation^. 

The next step will be to take the parametrized Minkowski action describing charged 
spinning particles and the electromagnetic field on arbitrary spacelike hypersurfaces, in which 
the minimal coupling describes correctly the magnetic couplings, and to add non-minimally 
these electric couplings of the spinning particles with the {now dynamical) electromagnetic 
field. Since we have not found the needed covariant modification of the Lagrangian, we 
limit ourselves to present the form of these couplings after the reduction to the Wigner 
hyperplanes, viz. in the rest-frame instant form of dynamics. 

Therefore, we will have a modified invariant mass for the isolated system plus the ordinary 
rest-frame conditions, namely the vanishing of the total 3-momentum of the isolated system 



5 See Ref. [12] and its rich bibliography for the special cases in which this is possible and for the 
lack of a mathematical justification of the validity of the Foldy-Wouthuysen transformation. 

6 It produces forms of the Dirac-like constraint in which there are only polynomials in the Grass- 
mann variables which after quantization become even operators. 

7 In particular the approximation will be valid as long as the electric fields are not too strong. 
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inside the Wigner hyperplane. This will produce a modification of the equations of motion. 
In particular the electromagnetic field equations will have a more complex particle source 
term, one dependent on the pseudoclassical spin variables. Therefore, the radiation gauge 
Lienard-Wiechert solutions of Ref. || will be modified. 

At this point we have only to repeat the whole procedure, identified in the previous 
paper on charged scalar particles, to get the semi-classical phase space version of the unique 
semi-classical Lienard-Wiechert potential and electric field. Then we will make the canonical 
reduction to only particle degrees of freedom by imposing conditions on the electromagnetic 
field in the radiation gauge such that it coincide with this Hamiltonian Lienard-Wiechert 
electromagnetic field. After the identification of the new particle canonical variables, the 
final form of the Hamiltonian will identify the relativistic form of the Darwin potential for 
spinning particles. 

In the 2-body case, the quantization of the lowest order in c~ 2 of the N=2 final Hamilto- 
nian has to be compared with the positive-energy sector of the 3-dimensional positive-energy 
reduction of the Bethe-Salpeter equation |14]-|18 for 2-body bound states with spin 1/2 con- 
stituents of arbitrary mass (hydrogen atom, positronium, muonium). We find that our Dar- 
win, spin-orbit and spin-spin terms reproduce exactly the standard Bethe-Salpeter result: 
this is the first time that, by making a relativistic separation of the center-of-mass motion 
and a subsequent study and canonical reduction of the relative motion, results coming from 
quantum field theory can be exactly reproduced. 

The paper is organized as follows. In Section II we study the pseudo-classical Foldy- 
Wouthuysen transformation applied to the manifestly Lorentz-covariant spinning particle in 
an external electromagnetic field. A review of the theory of the positive-energy spinning 
particle interacting with a dynamical electromagnetic field in the rest-frame instant form of 
dynamics is given in Section III. In Section IV we introduce the non-minimal coupling of 
the electric field to the positive-energy spinning particles on the Wigner hyperplane using 
the results of Section II. In Section V we give the equations of motion for the particles and 
the field in the radiation gauge on the Wigner hyperplane. In Section VI we determine the 
unique semi-classical Lienard-Wiechert solution, we find its phase space expression, we add 
second class constraints to eliminate the radiation field, we find the Dirac brackets and the 
new canonical basis for the particles. The semi-classical Hamiltonian with the Darwin and 
spin-dependent potentials for the A-body problem is found in Section VII. In Section VIII we 
study the 2-body problem: we find the semiclassical Hamiltonian for muonium-, hydrogen- 
and positronium-like systems and then we quantize its lowest order part, reproducing the 
Bethe-Salpeter result. In the Conclusions after some comments we delineate how the research 
program could be developed to treat the non-Abelian case of the quark model. 

In Appendix A there is the determination of some functions defined in Section II. In 
Appendix B, after a review on spacelike hypersurfaces, there is an attempt to extend the 
results of Section IV outside Wigner hyperplanes and outside the radiation gauge. In Ap- 
pendix C there is the computation of the electromagnetic energy and 3-momentum when 
the Lienard-Wiechert solution is inserted in them. In Appendix D there the evaluation of 
some of the potentials of Section VII. In Appendix E there is the summation of the 2-body 
rest energy to a closed form. 
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II. THE PSEUDO-CLASSICAL FOLDY-WOUTHUYSEN TRANSFORMATION 
OF THE MANIFESTLY LORENTZ CO VARIANT SPINNING PARTICLE. 



In Ref. JTT| there is the coupling of the manifestly Lorentz covariant spinning particle to 
external electromagnetic fields. It is based on the singular Lagrangian 
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which, besides the standard minimal coupling, has a non-minimal mass renormalization 
—ieF^^^ = eF^S^. The ie coefficient in front of F flu £ i tJi £ i u , corresponding to the absence 
of an anomalous magnetic moment of the electron, is the only one ensuring that the two 
constraints remain first class even in presence of an external electromagnetic field. 

Indeed, besides the second class constraints — « and the added one (like in the 
free case) 7r 5 + -£5 ~ |, one gets the first class constraints 

Xd = (p„ - eA^x))C - m£ 5 ^ 
X=(p~ eA(x)f - m 2 + ieF^ u (x)eC ~ 0, 



{Xd,XdY 



IX, 



{x,xV = {x,xd}* = 0. 



(2.2) 



As said in the Introduction, if we have many particles we assume that the spin Grassmann 
variables of each particle commute with those of the other particles. 

Following Ref. |19| , we now describe also the electric charge of each particle in a semi- 
classical way by means of a pair of complex conjugate Grassmann variables f] #j(r), (9*(t) 
l"9fl satisfying (J, = I* = 6*6i is the generator of the U em {l) group of particle i) 



ef = 6f = 0, 



m + eWi = 0, 



9i9j 



8j6i, 



a a* a*a 



6*6* 



6*6*, 



I ^ J, 



&i6*6i, 



(2.3) 



8 After the elimination of the Grassmann momenta by going to Dirac brackets with respect to these 
second class constraints, the original Poisson brackets {x^,p v } = {£^,n v } = —vf", {£5,^5} = — 1 
become the following non-null Dirac bracket: {x fl ,p u }* = —rf 11 , {£ >1 ,£ u }* = in^ u , {£5, £5}* = — i', 
in what follows we will use the notation {., .} for these Dirac brackets. 

9 They are assumed to commute with the spin Grassmann variables. 
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The action now depends also on the configuration variables 9i(r) and 9*(r),i = 
1,..,N, through an extra kinetic piece for the complex Grassmann charges J |[^*(r)^j(r) — 

e;(T)ft(r)]dT, 

The Grassmann momenta associated to these extra variables give rise to the second class 
constraints Hoi + |0* ~ 0, 7r e » 4 + |0j « 0; n ei and 7Tg. 4 are then eliminated with the help 
of Dirac brackets 

{A, B}* = {A, B} - i[{A, n ei + % -e*}{ve* i + -0,, 5} + {A, vr e , { + fofai + B}], 

(2.4) 

so that the remaining Grassmann variables have the fundamental Dirac brackets (which we 
will still denote {., .} for the sake of simplicity) 

{e i (r),e,(r)} = {e*(r),e*(r)} = 0, 

{9 i (r),9*(r)} = -i8 ij . (2.5) 



In Ref. [13 a pseudo-classical Foldy-Wouthuysen canonical transformation was intro- 



duced, which, after quantization, realizes an exact Foldy-Wouthuysen transformation of the 
Dirac equation in the case A (x) = 0, 9 gffi = 0. However in that paper the electric charge 
was not Grassmann-valued. 

r " is 



The infinitesimal generator S of the pseudo-classical FW transformation of Ref. |13 

S- = 2i(n.&60(a), (2.6) 
where [in the free case we get a = p 2 , 9(a) = 9(p)] 

P = I P I, 

u = p-qA(x), 

a = ti 2 -iQF hk (x)ee, 



V^(a) = (l/2)arctg(^j. (2.7) 



m 



The action T5 of the canonical transformation generated by S on a generic function F 
in phase space is 

T s : F - F + {F, S} + (1/2!){{F, S}, S} + .... (2.8) 



Let us now see what is the effect of the canonical transformation ( |2.8|) when the electric 
charge Q is Grassmann-valued, so that we have Q 2 = 0. If we consider its action on the 
Dirac-like constraint 

Xd = (po - QA {x))£ - (p- QA(x)) ■ i- m£ 5 « 0, (2.9) 
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{6>i(r),7r ei (r)} = {0?(r), ^(r)} = -<%; {vr ei + §0?,^ + 5%} = 
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we get from Ref. JT3| the result 



T s :(-IL-£- m£ 5 ) - -&Vm* + a. (2.10) 
We have now to evaluate the following two terms at order Q 

T s : -QA (x)Zo - -Q[A (x)£ + {A (x)£ o , S} + (1/2\){{A (x)£ o , S}, S} + ...., 

T s : Mo - Mo + {Mo, S} + (1/2!){{Mo, S}, S} + .... (2.11) 
Using £f = 0, we can write: 



{{...{QA (x)to,S},S},...},S} 



Q(2i) n $ n (p, x, f )&fo n - odd 



+ 0(Q 2 ), (2.12) 



Q(2i) n $ n (p, x,£)£o n-even 
with the odd functions Q n defined in Appendix A. Dropping the Q 2 term we have: 

T s : -QA (x)to -> -QFxfax&to - QGx(p,x,g)^ , (2.13) 
with the following two functions (they are odd, i.e. linear in the three Grassmann variables 

a 



Gl(p,X,£) = E^Lq (2fc+l + )| $ 2fc+l ip, x, . 



(2.14) 



Equally, we can write: 



{{...{Mo, S}, S},...}, S} 



Q(2i) n * n (p,x,OZ 5 Z 



n — odd 



(2.15) 



Q(2i) n ty n (p,x,£)£ n-even,n^0, 
with the odd functions \I/ n defined in Appendix A. Then we get 

T s ■ Po(x)£ -> Mo + QF 2 (p, x, £)£ + <3G 2 (p, x, |5f 5 f 0> 



(2.16) 



with the following two functions (they are odd, i.e. linear, in the three Grassmann variables 
I) 



In total, we can sum the two term to get 

T s : (po - QA (x))£ -> p + QF(p, x, f) £ + QG(p, x, |)f 5 fo, 



(2.17) 



(2.18) 
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where F = F 2 — F\ and G = G2 — G1. See Appendix A for the determination of the functions 
F and G. From Eqs. dA"26|) , fl&gp their final expression is [0{p) = (l/2p) arctg(p/m)] 

F(p,x,£) = A {x) -1 



(m + V"^ 2 + p 2 ) V m 2 + p* 2 
sin [2p0(p)] 



P 



L ^ F)r) u 



d9(p) dA (x) sin [2p0(p)] - 2p9(p) 



p ■ dA Q (x) p ■ £. 



(2.19) 



~f dp u dx u p 3 

The action of the canonical transformation Ts's on the Dirac-like constraint xd thus 
turns out to be 



T s : Xd -> Xd = Po + Q^(p, a;, £ - ^Va + m 2 + QG(p, 2, £)6£o- 



(2.20) 



The new constraint x'd contains some terms, linear in £0 or £5, corresponding to even diagonal 
terms in quantum case and some terms, linear in £5^0, corresponding to odd skew-diagonal 
terms, according to the rule: £ 5 — > 75,^0 — > 7s7o : 



X 5 75 



(po + QFYjo - Vm 2 + a + QG^jp 

skew— diagonal 



(2.21) 



diagonal 

In order to cancel the odd skew-diagonal terms at order Q, we define a unknown infinites- 
imal generator W = Qw(p,x,£) for a new canonical transformation Tw 

Tw ■ Xd Xd = Xd + Q{x'd> w} = 

= X'd + Q{Po£o - &yJm* + fP, w}. (2.22) 

By imposing the Ansatz: 

w = wi(p, x, 15^0 + w 2 (p> x, £)£ 5 , (2.23) 
the condition for the cancellation of the skew-diagonal terms is 



{Po£o - ^m 2 + P 2 ' w ^ x ' Ofo + w 2 (p, 06} + G(p, x, O^o = 0. (2.24) 
This condition is equivalent to a partial differential equation for the unknown functions 



w 1 ,w 2 : 



z)wiPo - (-i)w 2 Jm 2 + p 2 - {po, - { Jm? + p 2 ,Wi}£ 5 £ + G&£o = 



Wipo - w 2 y/m 2 +fp = 0, 

-{Po, w 2 } - {^/m 2 +p 2 , wi} + G = 0, 



u>2 



po 



■fa - . 1 p. 9u>! + G = 0, 



Po 9^1 



1 



\/ m 2 + p 2 <9x° v/p^+Tn 2 " 



p • 9wi + G = 0. 



(2.25) 



12 



By introducing the Green function g p (x) for p^d^ we find the following solution for wi, 
w 2 [see Appendix B for the tetrads e^ A {u{j>))} 

Po 



w 2 



y/ m 2 + p- 



■Wl, 



p^w^x,^ 



m 2 + fPG(p, x, £), 



4 



d yg P (x-y) dm 2 + p 2 G(p,?/,0> 



p»d»g p (x) = S 4 (x), 

0(±p^x fl )5 3 (e^(u(p))x ll ) 



x 



fp 2 



(2.26) 



This shows the existence of the canonical transformation Tw and the non-locality of the 
separation of positive and negative energies. Therefore, the odd term QGC, 5 C, of Eq.( [2.20D 
can be eliminated at the order Q with the canonical transformation Tw- 

After these two canonical transformations Tw o T$, we get the following form of the 
Dirac-like constraint at order Q 

1 



Xd 



p -QA Q {x)+iQ{p-i){i-E{x)) 



(m + y/m 2 + p 2 ) \/m 2 + p 2 



£ - Vm 2 + a£ 5 



[p - QA (x) - Q- 



p ■ E(x) x S 



£ - Vm 2 + a£ 5 « 0, 



^2.27) 



(m + y/ m 2 + p 2 )y/ m 2 + p 2 - 

where we introduced the spin S — — |£ x £. 

To find the final form x" {Tw °Ts ■ x ~~ > x") °f the mass-shell constraint we use the fact 
that Tw ° T$ is a canonical transformation, so that 

{x£,xS,} = iX ff - (2-28) 

We get (F hk = e hkr B r ) 

2 



X 



Po - QA (x) - Qp- E(x) x S 



m — \/m 2 + p 2 
p 2 yj m 2 + p 2 



m 2 + (p- QA(x)) 2 + 2QS ■ B(x) 



0. 



^2.29) 



after having discarded terms proportional to £o£s, since they vanish due to x'b ~ 0> which 
implies £5 weakly proportional to £ c at order Q. 

The mass-shell constraint can to be resolved at order Q in two constraints, corresponding 
to the two signs of the energy: 

1 



X 



(±) 



Po -QA (x)+iQ(p-0(t;-E(x)) 



(m + yjm 2 + p 2 ) \Jm 2 + p 2 



T \/m 2 + {p- QA{x)) 2 - iQti h ti k F hk {x) 

Qp-E(x) x S 

= p - QA (x) 



(m + \J m 2 + p 2 )y/ m 2 + p 2 



=F 



=F dm 2 + (p - QA(x)) 2 + 2QS ■ B(x) « 0. 



^2.30) 
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In this way we have identified the extra coupling at order Q to the electric field E{x) = 
— dgjo - 1 — 8Aq(x) of a spinning particle with a definite sign of the energy. Only the semi- 
classical approximation Q 2 = allows us to get this result in closed form. 

Now the problem is to interpret this result. Since Aq(x) and A(x) describe an external 
electromagnetic field, we are in an arbitrary fixed gauge. But we will need the result in 
the radiation gauge with both transverse radiation fields and act ion- at- a- distance Coulomb 
potentials from other charges acting simultaneously on the given charged spinning particle. 
It seems reasonable to interpret A (x) as the scalar Coulomb potential generated by the 

other charges. Then the term — ~ 9Aq(x) should be interpreted as the sum of the 
transverse radiation electric field and of the gradient of the Coulomb potential (action-at-a- 
distance electric field generated by the other charges), since this is the total electric field 
acting on the charged spinning particle. 
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III. THE POSITIVE-ENERGY SPINNING PARTICLE. 



In this Section we shall review the description of the positive-energy charged spinning 
particle given in Ref. |J. To define a charged spinning particle with a definite sign of the 
energy on spacelike hypersurfaces(^|, the starting point was the Lagrangian description of 
a charged scalar particle (see Refs. @J2(|) with only a real Grassmann 4- vector £ M (r) for 



the description of sphQ After the Legendre transformation to the Hamiltonian formalism, 
a Hamiltonian odd second class constraint is added to eliminate one of the components 
of In this way only three Grassmann variables will survive for each particle and, af- 
ter quantization, they will generate Pauli matrices acting on the 2-spinors describing the 
positive-energy wave functions. It turns out that the addition by hand of this constraint 
gives a consistent set of constraints. 

Since, due to this last constraint, the Lagrangian description is too complicated the 
model was defined directly in phase space by means of a set of constraints. As usual in 
relativistic particle mechanics, only the Hamiltonian description is tractable, because the 
Lagrangian one is too involved and very often it is impossible to get it in closed form. 

Given a 3+1 splitting of Minkowski spacetime with a foliation whose spacelike leaves are 
defined by the embeddings E i— > S r , (r, cr) i— > z^(t, a), the embeddings become new config- 
uration variables describing all possible hypersurfaces (all possible congruences of timelike 
accelerated observers). Since, due to the separate r- and a-reprametrization invariances 
of the action, there are first class constraints implying the independence of the descrip- 
tion from the choice of the 3+1 splitting, the z m (t, <t)'s are the gauge variables of this type 
of general covarianceL] See Appendix B for the definition of the notations and of the in- 
duced metric. Each positive energy particle is described by the coordinates f}i(r) such that 



ii 
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See Appendix B for some notions on spacelike hypersurfaces. 

To describe a positive-energy spinning particle we have to solve the first class constraint xd ~ 



of Eq.(1.3) to express £5 in terms of the Grassmann 4- vector The needed gauge fixing to xd ~ 
is a constraint eliminating one of the four components of £; M : as shown in Ref. || in the free case 
it is p^ + m£ 5 0, so that £5 ~ and p^ ~ hold simultaneously. The mass-shell constraint 
X ~ of Eq. Ql.l| ) is eliminated by the choice of the energy like for scalar particles. 

13 It is of the transversality type Pu£^ ~ 0, but with p^ being the conserved total momentum of 
the isolated system and not the particle momentum. 

14 In Eq.(29) of Ref. Q there is the Lagrangian generating all the constraints ( |3.3| ) except the 
transversality ones 0j(r) = p^£f ~ 0. The Lagrangian generating all the constraints fl3.3|) could 
be recovered by inverse Legendre transformation. This has been done in Ref. || [see its Eq.(50)] 
only in absence of electromagnetic field, because the general case is very complicated and not 
particularly interesting, except for the determination of the energy- momentum tensor. In any 
case, the Lagrangian contains both a minimal coupling of the positive-energy particles to the 
electromagnetic field and a non-minimal coupling like in Eq.(|2.1|). 
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The descriptions given by arbitrary congruences of timelike observers are gauge equivalent. 
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xf(r) = z^(t, fji^r)) . Moreover, for each particle there will be a Grassmann 4-vector £f(r) 
for the spin and a pair 9i(r), 6*{r) of complex scalar Grassmann variables for the electric 
charge. 

On the hypersurface E T , we describe the electromagnetic potential and field strength 
with Lorentz-scalar variables A a (t, a) and F a ^(t, a) respectively, which have the equal time 
concept introduced by the embedding, defined by [A, F^ are the standard electromagnetic 
potentials and filed strengths, respectively] 

A A (r, p) = z A (r *)Ap(z(T, a)), 

F M (r, a) = d A Az(r, a) - d s A A (r, a) = z£(t, a)z u s (r, a)F liV {z{r, a)). (3.1) 

The model of Ref. is defined in a phase space spanned by the canonically conjugate 
pairs of variables *"(r, a), p M (r, a) ; A A (r,a), 7r A (r,a) 16 ; ^(r), ^(r); £f(r), vrf(r); 0<(r), 
7Tej(r); 9*(t), vr »j(r) with the following Poisson brackets 

{z"{T,a),p v {r, a} = -r]^S 3 (a - a), 
{A A (r,a),^(T,a , )} = 45 3 (a-a'), 

{6»i(r),7r ej (r)} = -Sij, 
{9*(T),n e , J (r)} = -5 l3 , 

= -5^. (3.2) 

The total conserved 4-momentum of the system is = J d 3 ap^{r, a) [see Eqs.( |3.7| )]. 
The model is defined by the following set of constraints (only positive energy particles 
are considered) 

7r 9i + Z -9**0, 
ne*i + -9i ~ 0, 

Mr) = K(r) + l& (r)] / *o p„(r, a) = [vrf (r) + ^(r)] Pstl * 0, 
Xt(r) = <(r) - % -^{r) « 0, =► 0,,(r) » p^ftr) « 0, 

7r T (r, ct) « 

TV 

T(r, a) = Vfo <x) - Q t (r)5 3 (a - r^r)) « 0, 

i=l 

W M (r, a) = p Ai (r, a) - Z M (r, a) \ 1 7r r (r, a)gn(r, a)n s (r, a) + 



16 £'r = -Frr and £>r = 2 e fsi E si ( e fsl = eTSt ) are the electric and magnetic fields respectively; for 
9ab ~^ ^ab one S ets %r = ~ E ? = ET ■ 
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+ " l{ ] ,a) Y\r, a)j**(r, a)F fa (r, a)F s ,(r, a) + 



N 



i=i 



• y/rn$ - Y s (r, 3)(k, i¥ (t) - Q^A^t, 3))(k iS (t) - Q^A^r, 3)) 
1 N 

+ y , = E *V - ffi(r))S 3 (ff - %(r)) ■ 

Oi(T)Oi(r)^(T)^ (r)i 7 (r, ^iftjff (r)# (r)/«(r, a) 



+ 



+ 



<=i 

0*We?(r)tf(r) 



m 2 _ yjf^ ^^(^ACtf (t 
N 



l a (r, a)z Sf3 (T, a)7r s (r, «r) - 
Q<(r)ff(T)tf(r) 



- 2 E * 3 (- - m)-r-., ^— 

- • ^Jmf -7 rs (T, o-)k^(t)k w (t) 



»=i 



^ a (r, a)zvp(T, a)7 ru (r, ^^(t, B)F rs {r, 3) 



- Y s (t, a)z StM (r, a) [F^(t, a)vr"(r, a) + 

2V 

+ ^5 3 (ff-?/ ! (r))(K,H-Q ! (r)i f (r, a) 



0. 



(3.3) 



The Lagrangian density that generates all these constraints except for the combination 
4>i(r) = J d 3 a p^r, <r)£f (t) ~ of the third and fourth set above is given by [see Eq.(29) of 
Ref. HI 



N 



\mrMr) - ^(r)^(r)) - £Ur)#(r) 



m| - iQi(rmr)C!(r)z^(T, <f)z*(r, a)F M (r, 3) 
\J9tt(t, a) + 2g T ?(r, a)r]l(r) + g^{r, 3)rf: (r)?)|(r) 



Q i (r)(A r (r,5)+7)fA<(r,5))} 



^)%(^ ^O", <?)<?*>, a). (3.4) 



BO/ 



The following Dirac Hamiltonian (see Eqs.(46) and (47) of Ref. for its determination) 
completes the definition of our model 



H 



D 



d 6 a 



A"(r, a)K(r, a) - A T (r, a)T(r, 3) + /i T (r, 3)n r (r, 3) 
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N 

H;(t, a) = H,(t, a) + % £{W>> Xi(r)}xUr) + 

i=l 

i N 

~ ~ E{^(r, <?), Hr)}Ur) ~ 0. (3.5) 

Ps j=l 

See Ref. || for their Poisson brackets. The constraints 7r e j + |0* ps 0, 7r e »j + |^ ps 0, 
0,; ps 0, Xi ~ are second class, while the other are first class. 

The Dirac multipliers A M (r, a), /z T (r, (?) in the Dirac Hamiltonian imply that the con- 
straints 7i*(r, (?) ~ 0, 7r T (r, (?) ^ are primary constraints of the unknown Lagrangian; 
the other primary constraints are the second class ones, whose associated Dirac multipliers 
are determined by the Dirac algorithm as shown in Ref. ||. The Gauss law T(r, (?) ps is 
the only secondary constraints and it appears in the Dirac Hamiltonian with A t (t, a) as a 
multiplier not determined by the theory (it is a gauge variable). 

As said in Section II the Grassmann momenta ^(r), TTg*i(r) may be eliminated by using 
the Dirac brackets ( |2.5| ). 

The conserved Poincare generators are 



r N 
J"»= / d 3 a{z»(T,a)p»(T,a) - z»(T,a)p»(T,a)} - E [W<(r) - ^ (3.6) 
J i=i 

Since is a constant of the motion independent of the isolated system under investiga- 
tion, we may write 4>i — (^f + \^i)Ps^ ~ 0. Due to the transversality to the total conserved 
4-momentuinF], we have the possibility of reducing the £f s from 4 to 3 for each particle 
independently from the interactions. 

As we see, the component of 7^ M (r, (?) along P(r, (?) (i.e. orthogonal to E T ) contains the 
electromagnetic energy density and also spin-spin, spin- electric field and spin-magnetic field 
interactions^ All of them are necessary to get the first class property for these constraints. 
Instead the components of H^(r,a) along z?(t,(t) )(i.e. tangent to S r ) contain only the 
electromagnetic Poynting vector as with scalar particles H. 



A. The restriction to Wigner's hyperplanes 

As shown in Ref. , the restriction from arbitrary hypersurfaces S T to hyperplanes E# T 
is done by introducing the gauge- fixings by pairing up the first class constraints 7i*(r,a') 
with 



17 This is a highly non-local property, since the whole hypersurface S r is involved in the reduc- 
tion. Note that p 1 ^ weakly coincides with the total 4-momentum of the isolated system by using 
H„(T,a) psO. 

18 Only the last one will survive in the rest frame. 
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C(r, a) = z»(t, a) - xjf (r) - fiftr)*' « 0, 

{C(r,?),W;(r,^)} = -T l SF(ff-d'), (3.7) 
and the Dirac brackets which eliminate these variables are 
{A, B}* = {A,B} - J d 3 a{AC(r^)}{n;(r,a),B} + J d 3 a{A,H;(r,a)}{C(r,a),B}. 

(3.8) 

The hyperplane is now described by just 10 configuration variables: an origin x%(t) 
and the 6 independent degrees of freedom in an orthonormal tetrad &^(t) [b^rj^b 1 ^ = rj^] 
with b^ = l^, where / M is the r-independent normal to the hyperplane. We have Zf(r,a) = 
bf(r), z?{t, of) = x^(T)+b^(r)a f , gfs{r,a) = -6 fi , 7 w (r, a) = -5 fS , 7(t,<t) = detgn(T, of) = 
1. The non-vanishing Dirac brackets of the variables i/t, S^", A^, 7T 4 , £f , 7rJ are 

K(r),#r = -ir, 

W(r),^(r)}* = ^, 

{5r(r),5f(r)}* = C 7 7^f(r), 

{er(r),^(r)}* = -<^,. (3.9) 

While is the momentum conjugate to x%, the 6 independent momenta conjugate to 
the 6 degrees of freedom in the b^s are hidden in S^ u , which is a component of the angular 
momentum tensor 

= Lf + + Sf , 
5f = 6{f(r) | tfWV(r, a) - &£(r) | d 3 aa ¥ p^(r, a), 

N 

i=i 

{J'*", J Q/3 }* = C% a/3 .P S , {Lf , Lf }* = C^Lf, 

C>£# = + Tfirgrr ~ - (3. 10) 

Next, we eliminate the second class constraints Xi ( r ) ~ 0> <M r ) ~ with the new Dirac 
brackets 

{A, B}* D = {A, BY + i J2{A, xT(r)} V(XT M, £}* - ^ ER fc(r)}'{&(r), 5}*. (3.11) 

i=l Ps i=l 

Now we have 
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%(r)p aii = 0, 

V Ps 7 

AT 

5 f = -^erer, p.„st = o. (3.12) 
i=i 



However, now we get the following non-canonical Dirac brackets on E# T 



i=l 



Ps Ps 



P 



s 



{erw,ej(r)}L = <(»r - = <n^ y . (3.13) 

Ps 

In this way, we have eliminated the components of £f parallel to in a Lorentz-invariant 
way. The spin of each particle is described only by 3 Grassmann variables and the spin tensor 
S^ u satisfies a Weyssenhoff condition. The angular momentum tensor becomes 

jf» = + + Sf , 

JV 

^JJ*- V £tnva/3 jjyS p/iua/3 gt^S 

( t Ijlv ria/3-|* p/j,uaf3 Q'yS 

{Lf + 5f , Lf + sfy D = (Lf + sq 6 ) , 

p/J.ua/3 PsPs v a , PsPs fj, B _ PsPs v 3 _ PsPs u a (o -tA) 

/^s "s Ps 

Since by asking the time constancy of the gauge fixings fl3.7D we get fU X th {r,o) = 
X^{r) + X\{T)b u f {r)a f \ A>(r) = -a£(r), A^(r) = -A^(r) = § Ef - W \if), the Dirac 
Hamiltonian becomes 

H F D = ~X»(t)H,(t) - 1 -\^( t )H^{t) + | d 3 ( T[-A T (r,a)r(r,(T)+ ^(^5)^^,5)], (3.15) 

and we are left with only 12 first class constraints 
7r r (r, a) 0, 
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N 



r(r, a) = <9 r V(r, a) - £ 5 3 (<? - ^(r))Qi(r) « 0, 



i=i 



N 



i=i 



+ £Jm 2 + g^-Q^^T))]^ 



A? 



+ ;E*)-#)) 



0<WQi(r)CW^(r) 



m 2 + /t 2 (r) A /m 2 + /? 2 (r 



i=i ^/m? + /c?(r) 



A? 



-5E 



--■-1 Jmf + Kf(rJ 



-&«a(T)6«/3(r)F^(r,?7i(r)) j + 



+ M r ){ / dV[7? x S].(r, ct) + £ [/^(r) - Qi(r)A<(r, ^(r 

i=l 

W v [t) = &£(r) J d 3 aa r H"{T,a) - &£(r) / d z aa*W(T,a) = 



brjV8^j(r)^j(T)b Ta + 



1 N 

Af , 

+ £^V™f + feto " ft(r)A(r,«(T))] J + 
i=i 

+ i £ rf(r) gi( ; T) f (7 lf P bMr)Ar, Mr)) + 

i=l 



N 



- o E^( T ) / = = ha{T)hp{T)Fw{T,r]i{T))\ + 



+ (6jf(T)6Sf(T)-6?(r)6S?(r)) dW[* x 5]„(r, 5) + 



N 



+ E^(r)Mr) - Oi(r)X,(r,*(r))]} « 0. 



i=i 



(3.16) 



The next step [|[] is to select all the configurations of the isolated system which are 
timelike, namely with p 2 > 0. For them we can boost to the rest system with the standard 
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Wigner boost L^(p s ,p s ) for timelike Poincare orbits (see Appendix B) all the variables of 
the non-canonical basis x%(t), p^ , b^(r) , S pu (r) , A^(r,a) , n A {T,B), ffi(r) , /%(r) , 0i(r) , 
6*{t) , £f (r) with Lorentz indices (except p p ). This is a canonical transformation generated 
by e^'^^ with generating function 

Hr) = \u{Ps)IAPs)S» v {t), 



Kp) = II I{pYu II 
VPo 



"ft' 

P m ■ 



i 3 (p) = Hp), 



cosh u (p) 



sinhu;(p) = rj 



\P\ 



'p VP 
L^p,P) = ex P [u(p)I(p)t = 

= [ cosh (u(p)I(p)) + sinh {uj{p)I{p))Y v = 

= [1 - I 2 (p) + I 2 (p) coshcj(p) + I{p) sinhw(p)]^, 

Kihp) = exp [ - w(p)/(p)]^. 



(3.17) 



Since we have £fp S/1 = 0, we get I^ l/ (p s )S pu = 0, so that the addition of S£ u to S pu in JF is 
irrelevant. 

The new non-canonical basis (with the same Dirac brackets) is 



X s = X s - ~ e A U (Ps))VAB 



1 

— ( 

2 



i 



X 



M _ 



Vs^P 2 s(P° s + Vs^P 2 s ) 
1 



PS 



Vs\/P 2 s 



P°s+VsJp 2 VsJp 2 (P° s +VsJp 2 s) 



Ps Psi '!: 



k i — k 

e* = e*, e l = e 



A: 



7T 



S P<T = 



= sr + 



VsJp 2 (P° s + Vs\ P 2 S ) 



Ps^p: - rf) + v^WX - S^Ps)} , (3-18) 



where u p {p s ) = Ps/Vs\JP 2 — L%(p s ,p s ) [r] s = ±1; from now on we will select only the 
positive-energy branch ?y s = +1]. 

For later use, let us introduce the spin tensors 

Sf B = e A (u(ps))e B (u(p s ))S^, 



sr = ^(u(ps)K(^ps))s, 



22 



Since = implies £ ir = ^i^(p s ) = ^L^{P s ,p s ) = 0, we can reduce to 3 for each 

particle the Grassmann variables describing the spin 



We get 



£W := el(u(p s M(r) r = l,2,3. (3.20) 



Sf B = 0, 

AT AT 

i=i i=i 
i n 

Of _ ruv quv \ 1 or 

£ o 5 ^-^ »€' 

Z i=l 

= -^""W- (3.21) 



The £[( r )'s satisfy 



If we define 



W,G}h = - d€r i U{Pa)) %. (3-22) 



= \4{n{Vs)WB d -^^(S7 + S\% (3.23) 



we get 



{^,^}L = 0. (3.24) 

Therefore, with respect to the Dirac brackets {.,.}£> we have obtained a basis in which 
x^(t), p*£, Aa(t,<t), tt a (t,(j), ffi(r), /%(r), £[(t), 0i(r), 0*(t), are canonical variables and 
only b^(r), S^ u (r), are not canonical. The new canonical origin x% of the hyperplane has 
the same noncovariance of the Newton- Wigner position operator and describes the external 
decoupled canonical center of mass of the isolated system. In terms of this variable we get 
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,D£ JD — D £ ) 

J Q/J.V QCtf3\* q-yS 

l D s > D s ID — Ly 7 <5 °s J 

f qfiu q q /9\* r n/if fa/3T* f Ql Jt > 1 ' fa/3\* n 

i^s > D £ ID - Ps i^s SD — \ a i ) L s ID — U 5 

{J"", J a/3 }^ = C^P & . (3-25) 

As shown in Ref. ||, we can restrict ourselves to the Wigner hyperplanes £vkt ; orthogonal 
to p%, i.e. with normals P = u^{p s ), with the gauge-fixings 



rj(r) = 65(r) - 6^ =A ( M (p s )) « 

=► 6j(r) = <(«(p.))&S(t) « itf, (3.26) 

which imply the new Dirac brackets 

{A, BY* = {A,B}* D - ^[{A,H^r D ( Vja e^(u(p s )) - r )5a E°{u(p s )j){T° D ,B} D + 



+ {A,TB}h[Va^(u( Ps ))- VcT ^(u(p s ))){H^,Br D . (3.27) 



The time constancy of the gauge-fixings ( |3.26| ) implies A Mly (r) , &^(t) = L^a(Ps,P s ) 
and 7t^(t) = (namely the determination of S^" in terms of the variables of the system). 
The remaining variables form a canonical basis 

K(t)X}d = -rT, 

{A A (r, a), n B (r, a')}** = V *5 3 (a -a'). (3.28) 

As shown in Ref. Q], the dependence of the gauge-fixing ( |3.26| ) on p^ implies that the 
Lorentz-scalar indices A = {r,f} become Wigner indices A = {r, r}: i) A A=T (r,cx) is a 
Lorentz-scalar field; ii) AA= r (T, a), £[(t), ?7[(t) , are Wigner spin 1 3- vectors which 

transform with Wigner rotations under the action of external Minkowski Lorentz boosts. 

On Y*wt the Poincare generators are 

Fsi ° s ~ ^sPs ^sPs ' ) 

£0i = £ij = 5 ir 5 js S rs^ ^ 32g ^ 

because one can express S^ u in terms of S AB = e^(u(p s ))e B (u(p s ))S f11 ' . 
Since H^ v (r) = implies 
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Sf = (eJf(u(p.)K(p.) - e v MPs))u^p s )){\ J d 3 aa r (n 2 (r, a) + B 2 (r, a)) + 

+ f>[(r) vV 2 - iQi(TmT)$(T)F»v(T,fH(T)) + ft(r) - Q,l(r,^(r))] 2 } + 
i=i 

- (e?(u(p.))<(ii(p.)) - <(u(p.))e?(u(p.))){ J d 3 aa r (n x fl) a (r,<?) + 

N 

+ J2vKr)Mr)-Qi(T)A s (T,f] i ))}, (3.30) 



we get 



Sf B = (5*5* - 5*5*) [i / d 3 aa r (n 2 (r, a) + 5 2 (r, a)) + 



+ 5>rtoV™? - <Q < (T)ff(r)6'(T)F t(l ,(r,^(T)) + fo(r) - Q^(r, i*(t))]' 
i=i 

- (6?6f - 5?5f) [ J rfW(7? x B) s (r, a) + 
+ J2vXT)(k is (T)-Q i (T)A s (T,f) i (T))) 



+ 



(3.31) 



However, the boosts S° r = S' or [S'| r = 0] do not contribute to the previous realization of the 
Poincare generators: this is the external Poincare algebra in the rest-frame Wigner-covariant 
instant form of dynamics. 

The original variables z^(t,ct), p M (r, a), are reduced only to x^, p^ on the Wigner hy- 
perplane T> Wr . On it only 6 first class constraints survive 



7r T (r, a) pa 0, 
r(r, a) = <9 r ^(r, a) ~ E Q,,5 3 (a - ^(r)) » 0, 



iV 



i=l 



W(t) = - K(p s )^(r) + e^( M (p s ))H pr (r)] = 
= u»( Ps )H(t) + e^(p s ))ft pr (r) » 0, 



where 



W(r) = Jp 2 - U rd {r) =r] s Jp 2 s - 



N 



tt 2 + B 2 )(t, «?) + 



+ E " V™ 2 - iQiW^We/ (r)F rs (r : t£(t)) + [^(r) - Q t (r)A(r, ff^r))] 2 » 0, 



i=l 



71 X B 



N 

( r > s ) + E [M r ) - Qi(r)A r (r, r/i(r))j « 0, 



{W{r),Tr{r)Y* = J d 3 a{[u^p s )e^MPs))-u v (p s )e^u(p s ))y(r,a) + 
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- e?(u( Ps ))F rs (T,ay s (u(p s ))}T(T,a) « 0. (3.32) 

Let us remark that in TC(r) ~ the spin-spin and spin- electric field interactions have 
disappeared on T, Wt , i.e. in the inertial systems associated with the Wigner hyperplane. 
There is only the spin-magnetic field interaction 

- iQi^tFrsir, Ur)) = -2Q f ft* ■ B(r, fj t (r)), S l( = x £, (3.33) 

like in the non-relativistic Pauli equation. This term has been put inside the square roots 
of the kinetic terms in analogy with the results of Section II. 

Therefore, we get a kind of relativistic Pauli Hamiltonian describing the interaction of a 
positive-energy massive spinning particle belonging to the (|, 0) representation of the Lorentz 
group with the electromagnetic field, whose non-relativistic limit is the pseudo-classical form 
of the ordinary Pauli Hamiltonian. 

The constraints TC p (t) « identify the Wigner hyperplane Ejy T (whose embedding 
in Minkowski spacetime is z M (r, a) = x%(t) + e^{u{p s ))a r ) with the intrinsic rest frame 
(vanishing of the total Wigner spin 1 3-momentum of the isolated system) and mean that 
the conjugate 3-coordinate a = q + of the internal 3-center of mass of the isolated system on 
Svi/r is a gauge variable. The natural gauge-fixing for TC p (t) « is q + ~ 0; in this way the 
internal center of mass coincides with the origin of Hwt- x s( t ) = z ^i. T i & = 0)- 

On E^v the Dirac Hamiltonian becomes 

H D = \{t)H{t) - A(r) • H p (t) + Jd 3 cr[- A t (t, a)T(r, a) + /i t (t, a)vr r (r, a)] , (3.34) 

so that the external canonical 4-center of mass has a velocity parallel to namely it 
has no classical zitterbewegung just as with the Foldy-Wouthuysen mean position. 

See Ref. |21[] for the relativistic kinematics of the N-body problem in the rest-frame 
instant form of dynamics and Ref. [0 for the non-relativistic limit oi7i re i. 

B. Dirac's observables and equations of motion. 



As shown in Ref. |22| , the Dirac observables of the electromagnetic field are the transverse 
quantities A±_ r (r,a) , 7?±(t,<t), defined by the decomposition [A = — d 2 ] 

A r (r, a) = d r n em (r, a) + A ±r (r, a), d ■ A±(r, a) = 0, 

Qr N 

n r (r, a) = vr^r, a) + — [T(t, a) - £ Q^{r)5\a - ff^r)) 

t-^" i=i 



d 

Vem(r,a) = -— ■ A(r,a), (3.35) 

while the gauge variables are A T (r,a) and r) em (r,a), being conjugated to the first class 
constraints tt t (t, a) ~ 0, T(r, a) ~ 0. 

Concerning the particle variables, we have that «[(t), 9i(r), 0*(r), are not gauge invariant 
because 
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K(r), r(r, d)}% = Q^ 3 (a - jfc(r)), 

{^(r), T(r, a)}^ = ieA{r)5\B - ^(r)), 
{0*(r), r(r, = - #(t)). 



(3.36) 



However, the position variables ^(t) and the spin variables £f(r) are g au S e invariant. 
The Dirac observables for the particles are obtained through a dressing with a Coulomb 
cloud 



^(t) = e iei ^ m( - T ' ni) 9i(T), 
B*(t) = e - ie ^-( T >*)^(r), 

Ki(r) = Ki(r) - Qi{r)dr] em (T, a) k^t) - QiA±(r, ffi) = k^t) - QiA(r, 



The electric charges are gauge invariant 

Qi = e i e*(r)6 l (r) = Q l = e ?; ^(r)^(r) [4(r) = => Q<(r) = Q, 
When the Gauss law is satisfied, T(r,a) = 0, Eq.( |3.35| ) implies 

QiQj 



1...N 



d 3 an 2 , (t, a 



(3.37) 



(3.38) 



(3.39) 



so that we get 



H(r) 



h 2 s~ ftrel = 

-{lJd*o-(7ri(T,a)+B 2 (T,a 



QiQj 



1...N 

y — 

~^j ^ I Vi(T) - Vj(T) 



+ 



N 

+ E VW ~ iQi&(T)8(T)F r .(T, r£(r)) + fe(r) - Q t A ± (r, ^(r))] 2 } « 0. (3.40) 



i=l 



We see the emergence of the Coulomb potential from field theory and the regularization 
of the Coulomb self-energy (the rule) due to the Grassmann character of the electric 
charges, Qf = 0. In this way all the effects of order Q\ are eliminated, but not those of 



order QiQj, i ^ j |20[ . 

Let us remark that even if we have not fixed the electromagnetic gauge but simply de- 
coupled the gauge variables A T (r, a) and r/ em (r, a), this procedure is equivalent to a Wigner- 
covariant radiation gaugeQ. 

There is no odd first class constraint, because massive 2-spinors do not satisfy any spinor 
equation. The quantization of this Hamiltonian in the free case gives a non-local Schrodinger 
equation = ±V m ' 2 + ~^$±i with the kinetic square root operator [23], for a 2-spinor, 



19 The associated natural gauge fixings would be A t {t,3) and rj em (r,a) as 0. 
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which corresponds to the upper (or lower) part of positive (or negative) energy Dirac spinors 
boosted at rest 0. These 2-spinors are parity eigenstates in the rest frame. 
The three constraints defining the rest frame become 

N . r 

H p {t) = Y / K i (r)+ d 3 an ± (r, a) x B(r, a) « 0, (3.41) 
i=i 

and are independent of the interactions, as expected in an instant form of dynamics, like 
the 3-spin 



S7 = £[tf(rK(r) - (t)i$(t)] + J <Pa[o*(* ± * B) s - a*(n ± x B) r }(r, a), 

i=l 

N 

S rs = S r s s + S? = S r s s - i £ Cid- (3-42) 



i=i 

The Pauli-Lubanski 4- vector and the spin Poincare Casimir are [e s = 

yy s — g Psv^spa — 2 Psv^po — 

— * p . Cj 

= (Ps-S; e s S + — p s ) = 

= \^ pu Ps»{S sp « + h pa ) := W<P» + E?, 

= -\p 2 AJ» v = -P 2 S S. (3.43) 

-± -± o 

This shows that £ = S s + S% is the rest-frame Thomas spin: W Sfl (p s ) = W sy L y ^{j) s ]P s ) = 
(0;ej). 

As shown in Ref. Q, it is convenient to replace the external center-of-mass canonical 
coordinates x*£, p^ with a new basis defined by the following canonical transformation (T s is 
the Lorentz-scalar time of the rest frame) 



T = Ps tm e = J v 2 

J- S l , t s W p s 



K = %. (3.44) 



The inverse canonical transformation is 



20 So that they also coincide with the corresponding parts of the positive (or negative) energy 
Foldy-Wouthuysen spinors boosted at rest. 
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p° s = e.Vl + kl (3.45) 

By adding the gauge-fixing T s — r « f\ whose time constancy implies A(r) = — 1, 
we get a frozen phase space with the Dirac Hamiltonian Hp = — A(r) • 7^ p (r). Like in 
the Hamilton- Jacobi theory, to reintroduce an evolution in r = T s one uses the energy as 
HamiltonianQ so that the Dirac Hamiltonian for the rest-frame instant form of dynamics is 

H D = n rel (r)-X(T)-H p (r). (3.46) 

The effective Hamiltonian is the invariant mass M of the isolated system 

N 



M = H rel {r) = J2 VX 2 - mi(r)a(r)F rs (T, 7*(t)) + [£ t (r) - Q.A^r, r^r))] 2 + 



i=i 



1...N | 



N 



= E v m * ~ ' ^ + - 3^( r > ^)] 2 + 

i=l 

+ E m 4 , a(T) _ ffi(T) , + /^U^J (r , g) . (3.47) 



See Ref. [§] for the associated Hamilton equations. 

In the next Section we will modify this invariant mass to include the non-minimal cou- 
plings of the electric field identified in Section II. 



21 It identifies the rest-frame time T s with the parameter r of the foliation of Minkowski spacetime 
with the Wigner hyperplanes associated with the isolated system. 



22 



See Ref. [24] for a different demonstration of this result. 
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IV. NON-MINIMAL COUPLING TO THE ELECTRIC FIELD OF THE 
POSITIVE-ENERGY SPINNING PARTICLE. 



In this Section we modify the invariant mass ( 3.47|) on the Wigner hyperplane to include 
the non-minimal coupling to the electric field of Section II. Then we define the external and 
internal realizations of the Poincare group. 



A. The Non-Minimal Coupling to the Electric Field. 

As shown in the previous Section on the Wigner hyperplane of the rest-frame instant 
form we have the following four first class constraints for the description of positive-energy 
charged spinning particles coupled to the electromagnetic field in terms of Dirac observables 
corresponding to the Wigner- covariant radiation gauge 

WW =«. - {| fd^Kr.a) + B\r, 3 j) + f ^ , .ffi ... , + 



; 4tt I f}i(r) -%(r) 
+ E V™? - iQ&{r)8 (r)F rs (r, tfc(r)) + fe(r) - Q^r, ffi(r))f} « 0, 



i=l 
iV 



i=l 



W P (r) = £ *»(t) + / rfVvr ± (r, <?) x 5(r, a) « 0, (4.1) 



The original minimal and non-minimal couplings to the electromagnetic field produced 
the Coulomb potential and the non-minimal coupling to the magnetic field. Now we have to 
introduce the modification suggested by the pseudo-classical Foldy-Wouthuysen canonical 
transformation of Section II to include the semiclassical non-minimal coupling to the electric 
field. 

We see that the kinetic term of each particle is the same in Eq. (|4.1| ) and ( |2.30| ). In 



Eq.flOoD there was no choice of the electromagnetic gauge. Instead Eq. (|4.1| ) is in the 



radiation gauge, so that instead of the vector potential we now have the transverse vector 
potential of the radiation field. 

On each particle i acts: i) the radiation field; ii) the scalar Coulomb potential of the 
other charges. As said at the end of Section II, we have to replace Aq(x) of Eq. (|2.30|) with 
just the scalar Coulomb potential 

acting on particle i (the total Coulomb potential is V — hYliLiQiVi) while the term <f ■ 



[ ^ J) + 9v4o(x)) of Eq.( p.30|) will be replaced by two terms: i) the transverse electric field 



£ • 7?j_(t, o = ?7i(t)) of the radiation field at the position of the particle and ii) the spatial 
variation of the Coulomb field of the other particles when particle i moves appearing in the 
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combination^ 



^Q ^ ^(r)■e(r)e(r)■4^(r/ i (r)) 
(mi + yj m\ + ^Jr))\fmf + n t (r) 



(4.3) 



On the other hand, the rest-frame conditions 7Y p (r) m remain the same, since they do not 
depend on the interaction in an instant form of the dynamics. 

Therefore the modified constraints on the Wigner hyperplane, written in terms of the 
Dirac observables, are 



H = e s — M 



N 

e s -(EV m ?- iQiCI(r)a{r)Frs(r,ffi{r)) + *i(r) - Q^ ± (t, ^(r))) 2 

i=l 

»=i (mj + \/m 2 + k 4 2 (t)) \/m 2 + «?(r) 



r QiQj 



QiQjZiir) ■ &(t) &(t) • (^(r) - ^(r)) 



4tt I 7fc(r) - %(r) | 3 (m, + \/ m? + /?■ (r)) ^m 2 + kJ(t) 
+ / d 3 ai[^l + if](T,a 



+ 



at / : 



i=l 



. -QMr) x &(r) • B(t, ffi(T)) . Q^O) • &(r) &(r) ■ tt ± (t, ^(t)) 

+ 2 = = I 



z2 , 



2\/m 2 + £~(rl 

QiQj 



(rrii + yjmf + K^r^^mf + K t 2 (r) 



+ 



^ l 4tt I ^(t) - 7fe(T) | 

QiQjZiir) ■ &(r) &(r) ■ (^(t) - %(r)) 



4vr | f/i(r) - ^(r) | 3 (m, + \/ m? + k](t))\I mf + k-(t) 
+ / ^V-^l + S ](r,a)) «0, 



+ 



W P (r) = E^M+ / d 3 a7r ± (r, 3) x 5(r, 3) « 0. 

i=l 



(4.4) 



One can check that these constraints are still first class. 



23 As shown in Section VIII it is just the quantization of the consequences of this term that will 
produce the spin-orbit and Darwin terms in the final semiclassical Hamiltonian. 
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The Dirac Hamiltonian ( |3.46|) is replaced by 



H D = M' -\(t)-H p (t). (4.5) 

See Subsection 2 of Appendix B for an attempt to find the constraints before the restric- 
tion to the radiation gauge and on arbitrary spacelike hypersurfaces. Since we do not know 
the Lagrangian, we must use Hamiltonian methods 24 for determining the energy momentum 
tensor T AB and this requires the modified Dirac Hamiltonian on arbitrary hypersurfaces, 
which replaces Eq.(|3.5|), as input. 



B. External and internal Poincare Groups. 

The external realization of the Poincare algebra, given in Eq.( |3.29| ), is 



fj : .Kiri - + 5"'V s S rs , 



fitr crs„s 

Jf ■■ = *2p£ - - - — (4-6) 

Ps+Vs^JPj 

where S r ee S r s + 3$ = e™ ( E£i tf(r)*y(r) + / fa ^(tt ± x B) v (t, a) - f E^i &(t)?(t)) 
as a consequence of Eq. (|3.42j ) . 

The internal unfaithful realization of the Poincare algebra has the generators [from now 
on we denote by M the M of Eqs.( f4.4|) , (|4.5| ); the Green function c(cx) is defined in Eq. flB15|) 
of Appendix B] 



n / ; 



i=i 



.QM T ) x U T ) ■ B ( T , ^( r )) • Qi%i(T) ■ &(r) &(r) ■ n±(r, ?7i(r)) 

+ * , „ = I- 



2^/m 2 + K t 2 (r) (m, + \frnj+ k^t)) \j m 2 + £ 2 (r) 



EQtQj 



-4tt I rfi(r) - ^(r) 



24 If we have the action S = j dtL of a non-singular system with configuration variables q a (we 
consider a finite-dimensional case for the sake of simplicity) coupled to an external gravitational 
field, the energy momentum tensor is defined as T^ u (z) = — -j= s SS r z \ \g=rj- The canonical mo- 



/g(z) s 9^(z) 1 

menta and the Hamiltonian are p a = Jp^a and H = Jpf-g"— L, respectively. Note that the momenta 
depend on the gravitational field. Therefore, by using the notation Sh = J dtH = j dtp a q a — S, 
we get T^z) = -^=A^\ p \ g=r) - ' / dt[c? - = (-^A\p)\^ ^ 



using the first half of Hamilton equations. This gives the Hamiltonian form of the energy momen- 
tum tensor. In the case of singular systems, the Hamiltonian has to be replaced with the Dirac 
Hamiltonian: the Hamiltonian form of energy momentum tensor is now dependent explicitly on 
the Dirac multipliers. 
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— I 



— — * — * 

QiQjKi(r) ■ &(t) &(t) • (rfi(r) - ^(r)) 



+ 



Af / : 

^[^/ m 2 + ( ^. (T) _g.^ (r ^. (T)))2 _ 



j=l 



, . Qi^(r) x &(r) ■ B(T,fji(T)) _ . Qjkjjr) ■ &(r) &(r) ■ jf ± (r, ^(r)) 1 

"l - ^ I o ^ y o y o ~T~ 



2Vmf + ^(r) 



+ E[ 



rrii + ymf + /c. 



Utt I ffi(r) - ffjir) 
d 



(mi + \/ m? + ^(r)) V m f + «i ( r ) 
+ 



+ 



x&(t) • y dVi c(<ti - ^(r)) • c(<7i - ^(r)) 
/rf'Vi[^ + 5](r,a), 



AT 

P(int) = = E ^( T ) + / X ~ °> 

i=l J 



A 



^n t ) = ^ = ^ S ^ S * = EftOO >< ^(^)) r " o E& >< + 

i=l Z i=l 



+ 



1 dV(ffX [^xiffr,?), 



N 



K 



(int) 



oor or 



=1 



AT l.JV q 

+ E [ E ^^'[7^— C (^( r ) - ^(r)) - V r ,(r)cmr) - jfc(r))] + 



i=l 



+ Q^y rf 3 (T7rl(r, a)c{a- ^(t)) 



^2 



dV^ + B )(r,a) - 



V " x ^(r)-B(T,ff i (r)) .Q^F(t)^(t) ■ &(r) • #l(t, ffc(r))- 



i=l 



- *E 



2\[mY+ kJ(t) (m^ + ^/m? + /^(r))^/m? + /t-(r) 

QiQ^(r)^(T) • S(t) g(r) • (#(r) - ^(r)) 



»#j 4tt I r/i(r) - ^(r) | 3 (m, + ymf + kJ(t))V m? + k*(t) 



at 

= 5r = -5r=-E^)V m » ?+ ^( T )-^( T ' ^))] 2 + 
i=i 

AT l.JV I Q 

+ E [ E QiQil— «ZF^(r) - 77- (r)) - ^(rM^r) - ^(r))] + 



i=l 
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+ Qij d 3 afc r ± (r,a)c(a- fj^r)) 

N 



r / d 3 aa r (f L + B )(r,a) - 



+ 



i=i 



QiQjVi( T )^i( T ) -&( r ) 



(mi + \Jmj + K^T^mj + k- (r) 



Mi (m^ + + /^(r))^/m? + k*(t) 
-» dp 

d{r) ■ gz> J d 3 ai c(ai - ffi{r)) ■ c(ai - %(r)). 



(4.7) 



As said in Refs. |5||n|, the natural gauge fixings for the rest-frame conditions V(i nt ) ~ 
are Ku n t) ~ 0. These conditions identify the internal 3-center of energy of M0ller, which 
weakly coincides with the internal canonical 3-center of mass q + due to Vu n t) ~ 0. Therefore 
the gauge fixings JC^ int ) ~ are equivalent to q + ~ 0: they put the internal 3-center of mass 
in the origin x%(t) = z^ij, 0) of the 3-coordinates in each Wigner hyperplane. The gauge 
fixings q + ~ imply A(r) = 0. 
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V. THE EQUATIONS OF MOTION FOR N CHARGED POSITIVE-ENERGY 
SPINNING PARTICLES PLUS THE ELECTROMAGNETIC FIELD. 



The modified Dirac Hamiltonian ( |4.5| ), ( f4.7| ) in the rest-frame instant form on Wigner 
hyperplanes and with T s = r can be written in the form 



H d = M-\{t)-H p {t) = 

n / : 

= E [v m i + &(t) - q,a ± (t, for)))* + 



KAt) 



i=i 



2\ mj + K, t) 



QiQj 



1 - 



— i- 



+ I d 3 al[n ± + B}(T,a)-\(r)-n p (r) = 



+ 



A' 



E V m i + Mr) + fair) - QiA ± (r, ^(r))) 2 + 

QiQj 



i=i 



+ §4tt I ^(r)-^(r) 



— i- 



«i(r) • &(r) &(r) • (r/i(r) - ^-(r)) 



Vi(r) - %(r) | 2 (mi + Vmf + K*(r))^/m? + k^t) 
+ I d 3 ahf ± + B](r,a) - A(r) -7? p (r), 



+ 



(5.1) 



where we introduced the quantity 

Mr) = +iQS(r) x &(r) ■ £(t,?£(t)) - 2z 



• &(r) ■ 7fj.(r, 77,(7-)) 
rrii + \l mf + k, 2 (t) 



2Q, ^(r).5(r,^(r)) 



2Q i 5 iC (r). S(r,^(r)) + 



nx(T,fji(T)) x ^(r) 



^2, 



m, + ymf + /%(r) 
Analogously the modified Coulomb potential may be written as 



(5.2) 



y — 



1 + 



^( r ) - ^'( r ) I 2 + V m ? + ^( r ))\/ m ? + ( r ) 



(5.3) 
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By using the Hamilton-Dirac equations of motion, O = {0,Hc}, we obtain (where '=' 
means evaluated on the solutions of the equations of motion) 



Vi(t) = 



- A(r) - 



mf + TZt + (/?,(r) - QiA ± (t, r/,(r))) 2 

t j2 QiQMr)Ur)-(m(r)-vAr)) 

m 4tt I 77,(1-) - ^(r) | 3 (m, + yjmf + ^{r^y/mf + if^r) 

AT 



QiQA(t) -&(t)&(t) ■ (^(r) - 7fj(r))/c t (r) 
ft 4vr | ^(r) - rS-(r) |3 (m, + ^/m? + *?(r))(m t 2 + ^(r)) 3/2 
QiQjKijr) ■ fr(r) £i(t) ■ (rjijj) - ffjjr))^ = 

~« 4tt I 77, (r) - ?7j(r) | 3 (m, + \/m 2 + ^ (r)) 2 (m 2 + kJ(t)) 



+ 



+ 



-Q!i(Ml(r)) + l^ 



+ 



A(r) 



m, + (r) 



N 



QiQ^i(r)6(r)-(^(r)-^-(r)) 



j^i 47T | 77,(1-) - 77} (t) | 3 (m, + ymf + /^(r))\/m 2 + £, 2 (r) 



+ 



+ *£ 



Q/Q^(r) ■ fr(r) fr(r) • (g(r) - fyj^Wr) 
~i An | 77, (r) - 77} (r) | 3 (to, + yfm? + £ 2 (t))(to 2 + ^(r)) 3 / 2 

AT 



+ 



+ *£ 



QiQjKi(r) ■ &(t) &(t) ■ (ffi(T) - rjj{T))iZi{T) 



m 4vr | 77,(r) - t^t) | 3 (to, + ^ + ^(r)) 2 K + 



4 



and 



'to 2 + ^ (r) 
QiQk(Vi( r ) ~Vk(r)) 



Qi\ 1 = 



m 2 + re, (r) 



(5.4) 



+ 



m Air | ffi(r) - ff k (r) | 3 

-S(r)6(r) • f|^,(r) - ^-(r)| 2 / - 3(^(r) - ^(r))(^,(r) - ffj(r)) 



An I J/, (r) - ffjir) | 5 (to, + \jmj + K^r^mf + ^ (t) 
1 



+ 



Vm 2 + 7^(r) + (/?,(r) - Qzl±(r, ^(r))) 2 
(*i(r) - g^x(r,^W)«Q^^l(r,^(r)) + 
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y QiQk(Vi(r)-y k (r)) | QMr)u^A ± {T,fH(T)) | 
^ 4tt I jf, (r) - 77 fe (r) | 3 ^ + ^ (r) 



1 ggjjr] 

2 <9r? ; 



yJrnf + Ktir) 



.QA ■ 66 • [NOO - ^(r)| 2 / - 3(^(r) - ^-(r))(^(r) - ^-(r)) 



-Q*E 

i# 4tt I 77, (t) - ^-(r) | 5 (m, + V™/ 2 + ^i(^))\/ m l + ^( r ) 
To these Hamilton-Dirac equations we must add the rest-frame condition 

AT 



(5.5) 



E^) + / <P<t[t?± x S](r,5) « 0, 
i=i ^ 



(5.6) 



— * — * 

whose natural gauge fixing, implying A(r) = 0, is /C( int ) ~ 0, as mentionied in the previous 
Section. 

The equations of motion for the Grassmann variables and the spins are 



—i 



S(r) = 



j#i 4tt I if, (t) - 7fj(r) | 3 (m, + \Jmj + /? 2 (r))\/mf + (r) 

^(r) = Or(r)^(r), 

Jv ^ r JV 
W=EW = EW^(r), 



=1 



+ K. 



i=l 
jrst 



+ E 



e" t fl L (r,ifc(T)) + 



(mi + ^m, + R\)^mi + k?) 
(^(r)(vfi(r,^(r)) + i ^^(r)) 



^W^Kr^WJ + ^^^r)))]. (5.7) 



By using 



{i^(r, a), 7f ±s (r, ^(t))} = 5 rs 5 3 (a - $ (r)) , 



(5.8) 



the Hamilton-Dirac equations for the transverse fields are [P±((t) = S rs + Q-^r ls t ne trans- 
verse projector] 



TV 



A ±r (T, B) = - 7T±r(r, (?) - E ^WoO 



i=l 



iQi^(r)-e>)gW 3 (g-^(r)) 
(m^ + ymf + ^(r))Vmf + 



-[A(r)-9]i ±r (r )( r), 
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f r ± (r,a) = AA r ± (T,a) - [A(r) ■ d)ir r ± (T,a) 

■ QiH{r)Ur)-d8\a-^{r)) 



m, 2 + ^(t) 



-g<P['W(r)(y d (?-$(r)) 



^2^r AT 

- ^ iQiP± rs (v) [Ki(t) ■ ( 



«9r 2 



1=1 



(m^ + \/mf + K*(r))^/m? + At 2 (r) 



^ + Vmf + K-(r))(V / m J 2 + At 2 (r)) 3 



At,; T At,; T 



(m, + + ^i( r )) 2 (\/ m i + ^i( r )) 2 
(m^ + ymf + At 4 2 (r))\/^ + k?(t) 



^(r)^(r)5 3 (a-r/ 4 (r)) + 



+ 



^(r)-S(r)&(T)<P(*-i&(T)) 



(mi + Vmf + K?W)V m i + ^( r ) 
^i(r) • g(r) £(t)?,(t) • 45 3 (<? - # (r)) 



(mi + ymf + Ki(r))^/m? + At, 2 (r) 



(5.9) 



When we eliminate At/(r), ^(r) with Eqs.( |5.5| ), (|5.7|), by virtue of the Grassmann trun- 
cations due to Qf = we arrive at the following wave equation for the transverse potential 



d 2 A r ± (T,a) 
dr 2 



+ AA f ± (T,a)-[X(r)-d]7T r ± (T,a) 



— i- 



J2Qi[P r ±(m(r)S 3 (v-Vi(T))- 

i 

f ff(r)g(r) ■ &5\a - jfc(r)) 



+ 



+ P?(B)i 



\Jm 2 + At 2 (r) 
Ur) ■ ii(r) £f (t)^(t) • 4<* 3 (<? - tj& (r)) 



(mi + ymf + At 2 (r)) ymf + £ 2 (r) 
Expressed in terms of particle velocities the right hand side becomes 



(5.10) 



/(l-fe(r)g(r)-^ 3 (g-ffe(r)) 

m, ; 



1 - ^ (r) ^(r) • 6(t)^ (r)^(r) ■ <9 CT 5 3 (a - $ (r)) 



mi(^l-Vi(T) + 1) 



Let us now choose a gauge (the natural one is K.u nt \ 
gauge we obtain 



(5.11) 



0) implying A(r) = 0. In this 
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□Ai(T,<7) = (|^-^Ml(r,a 



dr 2 

N 



^(^■^^(r)-^(r)4 s (r)Vl-^(r) 



■3, 



i=l 



(1 + Vl-^ (r))m 



^r)Ur).d^l- Vl (r) + , + _ 
rrii J A 



TV 



£(t, a) = P['(<?) £ g,V|(r)5 3 (a - #(t)). 



(5.12) 



i=l 



We see that the effective current j±(r, a) has the structure of a pole-dipole. It is convenient 
to define an effective particle velocity operator (using boldface) V, 



Vf(r) = V°(t) 



mi 



1^, 



m;(l + V 1 — 77^ (r)) 



4 



QiVi(r) 



mi + Ks (t) 



m,- 



+ 



|^(r)| gjg(r) x £j(r) ^(r) ■ d CT 



:\Jm 2 + ^ 2 (r)(|^(r)| + ^/mf + ^(r)) 



(5.13) 



In the last line the operators Vj have been reexpressed in terms of the momenta by using 
Eqs.fl5ip. 

We point out that defining &(r) = r/,(r) = ^j&2 = and 4?° = d h ^jdr h and 

writing the particle equations of motion as (no sum over i) 



d A(r) 

-rim , — = 

dr . h an 



rrii 



QiF{ 



we obtain 



m,- 



(1-A 2 (r)) 3 / 2 



so that 



(5.14) 



(5.15) 



Qi(Fi-Mi'Fi)- 



(5.16) 



25 r = ci, %(i) = Vii T )> though using everywhere c = 1 , we have momentarily re-introduced it. 



39 



Thus in general we will have for every h > 1 

/3j QiG{. 



(5.17) 



so that using the Grassmann property of the charges 



o. 



h> 1. 



(5.18) 



This will lead to important simplifications later, allowing us to drop acceleration depen- 
dent terms in the force. 

Due to the projector P™(a) required by the rest-frame radiation gauge, the sources of 
the transverse (Wigner spin 1) vector potential becomes non-local and one has a system of 
integrodifferential equations (like the equations generated by Fokker- Tetrode-type actions) 
for which it is not known how to define an initial value problem. 

Let us end this Section with some comments on the equations of motion of the spin. 
Eq.( |5.12| ) shows that, besides the standard term for scalar particles |2Ti ], the particle current 

5 3 (a- 



contains also a dipole term P™(a) J2iLi (- 



771? +K; 



S , i4 (r)XK i (r) k;(t)-<9 



ffi(r)) in accord with the fact that the spinning particle has a pole-dipole structure [25 
according to Papapetrou's classification (see Refs. [27] for other pole-dipole models and 
Ref. [28] for their influence on the energy momentum tensor of action-at-a-distance models). 
For N=l (f/j i — > Tj^ vrti i — > m,...) we have 0: 

i) S? = e rst St = -iff', S r s s « Jd 3 a[((T r - rf )(vr ± x B) s - (o s - r/ s )(vf ± x 5) r ](r,a); 



ii 



-TSt Ct 

: ^ 



, 2 



The Bargmann-Michel-Telegdi equation |29[ for the spin part of the Pauli-Lubanski 4- 
vector of Eq.(|3^3|) 



V — TP 

m 



(5.19) 



for a spinning particle in an external electromagnetic field is replaced by the following 
equation in the canonical realization ( |3.29| ) of the external Poincare group in the rest- 



frame instant form for the isolated system of a positive-energy spinning particle plus the 
electromagnetic field (ft^ = 0) 



sir s 



1 



§ ir p s 



-f S ir S JS -e ■ 

2 ^ ^ uol poj rt 



(5.20) 



In the last line we used Eqs. (|5.7p for the spin. 



26 See Eq.(^43|) for the definition of 



2 e IJ-vpa 



P v SP a with PP^O. 
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VI. THE SEMI-CLASSICAL LIENARD-WIECHERT SOLUTION, THE SECOND 
CLASS CONSTRAINTS AND THE REDUCED PHASE SPACE. 



In this Section we study the retarded, advanced and symmetric Lienard-Wiechert solu- 
tions of Eqs.( |5.12|) . After an equal time development of the delay, we show that, by using the 
particle equations of motion, the higher order accelerations disappear by virtue of Q 2 = 0. 
Therefore, at the semi-classical level there is a unique Lienard-Wiechert solution. Then we 
obtain the phase space expression of the semi-classical solution. This allows to eliminate 
the radiation field degrees of freedom by adding second class constraints and to arrive at a 
reduced phase space containing only particle degrees of freedom. The result is a canonical 
basis of this reduced phase space. 

To simplify the notation from now on we shall denote the Dirac observables A±,... as 
A 



A. The Lienard-Wiechert Solutions, the Equal Time Expansion of the Delay and the 

Semi-Classical Approximation. 

The symmetric solution of Eq. ( |5.12| ) is 

A r ± s(r,a)= 1 -[A* ±+ + A* ± _](r,a) = 
I N O r 

= 2 E J d ^ d ^ P lV) V*(n)[0(r - n) + 9{t - r x )\ 

i=l 

S[(t - n) 2 -(3- 3 1 ) 2 ]5 3 (3 1 - ffiin)) = 

N O r 1 
= E ^ J dhPrWVKnWt - h) 2 --(a- fj^ch)) 2 ] : = 

N 

■ = T,Q^ Ar ±s ^ (r,a), (6.1) 
i=i 



in which we have put r = ct, P%{r) = ff^r) = - dv ^ and A± + = A±ret (Al- = AladvO 
for the retarded (advanced) solution. The equation for ti is c 2 (t — ti) 2 = (3 — ^(cti)) 2 with 
the two solutions being 

11 T 

t i+ {T, 3) = -T i+ {T, 3) = t - -r i+ (r l+ (r, 3), 3) = T l+ (r, a), 

c c c 

11 r 

ti-(r, 3) = -Ti-(r, 3) = t + -^--(^-(r, 3), 3) = - + 2i_(r, a), (6.2) 
c c c 

for the retarded and for the advanced case respectively. The light cone delta function is 

5[(r - n ) 2 -(3- ff,^)) 2 } = - h) 2 - I(* - ff^)) 2 } = 

5[ri_ - T i+ (r, 3)} + a[n_ -r,_(r,a)] ^ 



2|r-n -A(n) • (3-ff i (r 1 ))\ 2|r-n- ■ (5-^(n))| 
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The relative space location between the field point and the retarded or advanced particle 
position is 

<? - Vi{n±{r, a)) = fi±(r i± (r, B),a)= r i± (r i± (r, a), a)r i± (T i± (T, a), a), (6.4) 
and its length is related to the time interval by 

ri±{Ti±(r,a),a) = \ a - ^(r»±(r, a))\ = cT i± (r,a) = |r - r i± (r, 

=> r - r i± (r, a) = ±cT i± (r, a) = ±r i± (r i± (r, a), a). (6.5) 

The effective spatial interval is defined by 

Pi±(r i± (r,a),a) = ri ± (r i± (r, a), a)[l =p A(r i± (r, a)) ■ f i± (r i± (r, a), a)). (6.6) 

In terms of these variables, the retarded, advanced and time symmetric solutions are 

N O 1 
A r ±± (r, a) = Pr(a) £ ^ V*(r i± (r, ff)) 



i=l 



N 



i=l 



* 



1 



i=l 



j7T 



p i+ (r i+ (r, a), a) 



+ V|(r i _(r,a)) 



Pi-(ri-(r, #),&)■ 



We use the Smart- Wintner expansion [30 



f(n±) = f(r =F cT i± (T i± ((T, a), a))) = f(r - [±r i± (r i± (r, a), a))) 

oo I ^ d k-l 



fir) + E 



oo ( d k 



df(r 
dr 



(±r,fo*))' 



E 

fc=0 



A;! <ir fe 



f(r) (±r i (r,a)r[lTPi(r)-h(r,B)) 



where 



fi(r, a) = n(r, a)h(T, a) = a - ffi(r) = f i± (r ± (r, a), a)\ n ±(T,a)-- 



(6.7) 



(6.8) 



Pi±[Ti±) 



i±(Ti±)[l =F A(ri±) ■ r i± (r, a)] 



(6.9) 



The last line in Eq. (|6.8| ) is identical to the previous one since dr '^' cr ' > = — ' ^i( r ^). 
Hence we get 



Ai ± (r,a) = iT(^)E^E 1 |r^[^(^)^~ 1 (^^)] ) a\ s --{a a- ± _). (o.io) 
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In order to evaluate the above derivatives we need the Leibnitz formula for the kih 
derivative of the product f(r)g(r) 

d k . . f k\ d m f d k - m g 

dr kU9) J^ Q m\{k-m)\dT m dr k ~ m1 1 } 



Thus 



, n k\ dr kL 1 1 1 ^ fc! ^ m!(£;-ra)! dr fc - m dr m 

oo oo ('-r^ jk~m-\TS Am k— 1 oo oo / \/i+m Ah\rs Am h+m—1 

y^ y^ v+v a v j" 'i _ y^ y^ _L_ _j (q \2) 

^t m m\{k-m)\ dr k ~™ dr™ ~ m\h\ dr h dr™ 



Using the notation V^ 5 = we obtain the following expression for the vector 

potential 



Al ± (r, a) = P]_ s {a) Et^E ^Vf» s (r)0, ± ,,(r, a), (6.13) 

i=l 47F h=0 n - 



in which 

00 (m) m rfmji+rn-l/ -t\ oo /^Nm Jn„ _ . 

= E ^gr ^^ ( ' } = E ^^[^-^W) 2 ]^- 1 - (6.14) 



m=0 ml dr m ^ Q ml dr m 

In order to display the result of the evaluation of the derivative we use the formula 

d^ R{f{r)) - EE ni!naL ^ ( 2!^^ } - (6 - 15) 

(with the summations restricted so that X)r n r = ^, J2 r rn r — m ) to obtain 

00 ( T ) m ™ m! 



(/>i±,h(r,a) = E — r E E 



ml ^nf^. n ^ 1 



drj 1! I 2! 



(6.16) 



In this expression the symbol o represents a scalar product between the tensors to the 
left and to the right with the summation J2r n r = n indicating how the indices would be 
matched. Changing the m summation index to k = m — n we obtain 

00 00 / \k+n( \Ti r n r =n 

o,:,ir.o) EE E ( + ) 1 ) 



ra=0fc=0«i«2.. 71^2- ■■ 

d n r k+n+h -\r,a) ( &{t)Y 1 ( '/fV)^ 



o 



dfj VI!/ V 2! 



(6.17) 



(in this latter summation J2 r n r — n , J2 r rn r —n + k.) 
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Now we can take advantage of the Grassmann charges to significantly simplify the above 
multi-summations. As we have seen above, with a semi-classical Qi there are no accelerations 
on shell [Qify h ' = from Eq. (|5.18 )1 in the equations of motion of the particle since both 



the Coulomb potential and the Lienard-Wiechert Lorentz force on particle 'z' produced by the 
other particles, i.e. Qi[E±(r, fji{r)) + A( r ) x B(r, ffi(r))), are proportional to Qi. Therefore, 
the full set of Hamilton equations ( |5.5|) , ( |5.6|) for both fields and particles imply that at the 
semi-classical level we have a natural order reduction of the final particle equation of motion 
in the Lienard-Wiechert sector (only second order differential equations). One effect of this 
truncation is the elimination of multi-particle forces; all the interactions will be pairwise, 
in both the Lagrangian and Hamiltonian formalisms. This was to be expected since the 
rest-frame instant form is an equal-time description of the N particle system: (acceleration- 
independent) 3-body,.. N-body forces appear as soon as we go to a description with no 



concept of equal time, like in the standard approach with iV first class constraints |44] . 

Thus the only contributing indices are U2 = % = .. = 0, n\ = n and our expression for 
the transverse vector potentials simplify to 

N O 

A r ±± (r, a) = Pl s {a) £ ^V*(t)^ ± , (t, 3) = 

±[ A ] h n\ drj 

Al s (r,a) = \(Al + + A^_)(r,a). (6.18) 

'^2 



Since r« = Jff , we see that for odd n = 2m + 1 we get 

fl2m+l a2m+l 

' % 2m = -sgsmtfr = 0, (6.19) 




and this implies the equality of the retarded, advanced and symmetric Lienard-Wiechert 
potentials on-shell 

N O 

Al s (T,a)=Al ± (r,a)= P?(a) £ g^(r) 

^7ife4rt4 (6-20) 



% (2m)! V drf 



Therefore, at the semi-classical level there is only one Lienard-Wiechert sector with a 
uniquely determined standard action-at-a-distance interaction. 

B. The Phase Space Expression of the Semi-Classical Lienard-Wiechert Solution. 

We use a tensor notation to write the transverse symmetric vector potential above as 
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Using the definition of the Coulomb projection operator 



Px_ hk {B)F{a) = 5 hk F(a) - i- J d 3 a' ^ 



da h da k \a' - a 



■F(a') 



and compactifying the notation still further we obtain 



N Qi ^ l 



^i{r){m-daY m )\^-m\ 



\2m-l 



47T J \(T — (7 



\2m-l 



(6.22) 



(6.23) 



Integration by parts and changing from -J^ to -§g and translation gives 



N /n oo 



Qi ^ l 



-4vr^ (2m)! 



A ± 5(r, ?)=E?Et^ V 4 (r)(^(r) • <9 CT ) 2m ) |a - #(r)| 



2m-l 



,2m 



(a' - (a-^(r))| 



/2m-l 



a 



(6.24) 



We remind ourselves that, in addition to a convective part, the generalized current include 
magnetic and electric parts 



y.( r ) = ^.( r ) _ !Vj -Vi(r)Ct(r%(r) • d a + V 1 - ^?(r)^(r) ■ &(r)ff (r^Qr) ■ <j 



mi 



m;(l + VI -^i(r)) 



in which the derivatives act on 5. Due to the delta function in Eq. ( |5.12j) , this derivative can 
be replaced with minus the derivative on ffi(r). 

The integral in Eq. ( |6.24j ) is finite, and thus we can view it as the A — ► oo limit of an 
integral with a cutoff A and take the derivatives out. The integral is thus of the form 



-a CT (V,(r) • d a M(r) ■ d a Y m I rfV— 



a 



/2m- 1 



and 



— / d 3 -' 



a 



,2m-l 



a 



An J a \a'—(a — ffi)\ 2 Jo 



1 r A , / /2m+l 

da a 



\a' - {a-rfi{T))\ 



dz 



(6.26) 



1 f A , i ,2m+l 

da a 



2 Jo 



1 - 1 a' + {a — fji) 2 — 2a'\a — fji\z 
1 



a'\a - rfi 



a' + (a - rfi) 2 - 2a\a - ffi 



2\a — fj 

A 2m+l 



— — / da'a' 2m (\a' — \a — ffA \ — \a + \a 
-riAJo 



\a - 77, 



\2m+l 



2m + 1 (2m + l)(2m + 2)' 
Note that the A cutoff will get killed by the a derivatives. Thus, we obtain 



(6.27) 
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1 



2m 



\2m-l 



(2m + 2)! 



d a {Y t {r)-d a m{r)-d a Y m \a-^{r 



\2m+l 



N 



■ = E QiA±Si(Z - ■= A ±1 (a, r) + A ±2 (a, r) + A ±3 (a, r). (6.28) 



i=i 



We have introduced the notations 



N Qi ^ 1 



\2m-l 



A±i(r, ?)=EJE ^ V 4 (r)(^(r) • 9 CT ) 2m |<? - ^(r 
7 , - / \ f Q» r n | ■ g(r)g(r) • & \jl - %(r) " r 1 

AL2(T, (7) = (-) ^ — (1 + ^ 

*=i 47r mi(l + y/l -fh(T)) m =o L (2m + 2)! 



,2m+l 



\2m+l 



(6.29) 



and 



A± 3 (r, tr) = - E — ~ E h 

m= 



7i 4tt 



m, ; 



^ L (2m + 2)! 



9 CT (^(r)-9 CT ) 



^ T ).d a y m \a-^T 



|2m+li 



Note that since 



Ur)-d (T )Ur)-d t )F(\a-ff t (r)\) = 0, 
the A_|_ 3 (r, <?) term is zero. 



(6.30) 
(6.31) 



Using the result that (ff t • d,) 2 "^ 2 ™- 1 = [(2m - 1)!!] 2 ^[^ - (% • ^) 2 ] m [ 3, = 
we get 



r , ^ Aft y r [(2m -!)!!]• j , 



- 1 Mil 2 



°" - ^ n2 



I " - »7i 



yy 



(6.32) 



We point out that all three currents contribute to A±i, while only the convective and 
spin-electric parts contribute to Aj_ 2 . Using 

[(2m - l)!!] 2 (2m)! (m - 1/2)! _ v ^F(-) m - 1 (m - 1/2) 



(2m)! (m!) 2 2 2m 



(1/2 -m)\m\ 



v^(-y 



(-1/2 -m)\m\ 



= (-)* 



-1/2^ 



m 



(6.33) 



we find that 



AT r\ i oo 

^i(r,-) = E^v^)^ i ^E(-) 



N O ■ - 



l^-^l m=0 
1 



- (« ■ £4/ 



m 



I " - »7« 



i=i 



o- -77i 



(6.34) 
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In terms of canonically conjugate variables this becomes 



A±i[r, (t) = V -r- «i - «?i(r)§i(r) • <9 a + 

47T 



i=l 



7717 + I 771,; + \/m,- + ft 



1 



(6.35) 



For A^2(t, o)i we see from Eq. (|6.25 ) that we need an expression for (fj^r) ■ d v ) 2m+l \7]\ 2m+l 
and £i(r) • d^ff^r) ■ d ri ) 2rn \fj\ 2rn+1 . One can show by an induction procedure that 



| 2m+1 = [(2m+l)!!] 2 5:(-)^7j(^)) r 



\S-ih\' 



21 + 1 



(6.36) 



so that we get 



A± 2 (r, <t) = -) 2^ ^— C 1 / K 



•2 



i=l 



mi(l + Vl-^(r)) 



[(2m + 1)!!] 2 (-J 



£o (2m + 2)! 



Z=0 



r 



2/ + 1 



(6.37) 



Since we have 



[(2m+l)!!] 2 _ +1 /-l/2^ 
(2m + 2)! 1 ' lm+1, 



(6.38) 



we obtain (where I is the unit dyad) 



N n 

) = -V^*(i 



i=l 



mf + nf 



- 1 



mf + R 2 (mi + Jmf + Rf) \ a % 



m? + k 2 



K 2 - (Hi- p%) 2) 



a -rjj a -rjj 
\a - f)i\ \a - ffi 



(6.39) 



After combining some terms in Aj_ 2 (t, a), our exact result for the vector potential is 



N Q 

i=l 



mf + K-(m 8 + \/mf + Rf 



(1 



m 2 + Rf(mi + Jmf + R 2 ) \ a Vi\ 



a - 7]i a-7]i 



mf + k 2 



ttll* - ^mf + ft ■ ^ym 2 + Rf + yV 2 + ft ■ |§E^) 2 ) 
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O -Tfi O Tj. 



1 



iKi(r) ■ d a Ki(r) ■ &(t)&(t) 



■Ki ■ [1 — — 



(mi + yjm% + k?) 



(6.40) 



The corresponding E±s(t, a) and B±s(t, a) fields have the series forms 

i=l ^ n m=0 



1 



(2m + 2)! 



<9 CT (U,(r) • ^)(^(r) • 9 CT ) 2m+1 \a - Ur)\ 2m+l 



(6.41) 



where 



Ui(T) = ^(T) 



i\fl-Vi(T)Vi(T) ■ ii(T)ii(r) ffiir) ■ d a 



m 



(i + y/i-vdr)) 



(6.42) 



and 



A' 



2m- 1 



S 5 (r, c?) = - E S E 7^T V ^) x d M(r) ■ 3 a f m \a - ff^r) 

i=l 47r m=0 l Zm J- 

The corresponding esact forms in terms of the canonical variables are 

#±s(r , (Jj - - 2^ T- Kj - • Oa+ T== =; / 5 ^, 

j=i 4?r L ^/m 2 + ^(mj + y'm 2 + £ 2 ) 

r 1 1 



(6.43) 



(6.44) 



and 



7Tj_s(t, 5) = £_l<?(t, 5) = - E 



7 Qj 



1)- 



<7-77j| 



(6.45) 
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For the evaluation of the field energy we need Ej_ s + Bg, while for the field 3-momentum 
we need E±$ x Bg. Since the exact forms cannot be integrated to closed forms, we need the 
series forms for the fields. From Appendix D we get [^(t) = %-(T)?7y(r) = rji{r) — %(r); 
dij = d/dffij] 



and 



1..N 



U(r) := UdMEis + B 2 s )(T,ff) := E 

^ QiQj fs fs m ^ (^■4-) 2m+1 (^-4-) 2n+1 < n+2m+1 

k 4vr ^ ^J 1 3 (2n + 2m + 2)! 

rO. . A .VTT . . ft. .\( il. . ft. .\2m+l(Jf, . a. \2n+l„2n+2m+3 



+ 



+ 



(2n + 2m + 4)! 


— : + 


(Ui-4)(Vi-4)(^-4) 2m+1 (^ 


S \2n+l„2n+2m+3 


(2n + 2m + 4)! 






S \2n+l 2n+2m+3 
' u ij) '/ 



(2n + 2m + 4)! 
+ V- • V v " ' J ' v 11 u ij'h 



' r 3 (2n + 2m + 2)! 



(V, • ^)(V, • • fy)*"^ • ^) 2 ^g 



2n^,2n+2m+l 



(2n + 2m + 2)! 



/ d 3 a(E ±s x B s )(r,a) := ^ QiQjhi(rfi, rfj, Vij 
J i<j 



- v v (B v v ^-^) 2m+1 ^-^) 2n ^ +2 " ;+1 

4vr V «l * ' (2n + 2m + 2)! 



(u, • <%)(v, • • <%) 2m+1 (^ • d %j ) 2n n 



2n„ n 2n+2m+3 



(2n + 2m + 4)! 

Vj(y< - Uj) ■ %g ■ 4-) 2m+1 (j- • 4-) 2ra <" +2 " + 1 



(6.46) 



(2n + 2m + 2)! 

(i<— (6.47) 



C. Final Dirac brackets and their Darboux basis 



Till now we have worked in the reduced phase space of N positive-energy charged spinning 
particles plus the transverse electromagnetic field. This is a well defined isolated system with 
a global Darboux basis [ fji, Ki, A±(t, a), 7?j_(t, a)] and a well defined physical Hamiltonian, 
the invariant mass M = VJ int \ of Eqs. (|4.7| ), (|5.1| ). All possible configurations of motion take 
place in this reduced phase space. The space of solutions of Hamilton's equations is a 
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symplectic space, since there is a definition of Poisson brackets on the space of solutions. In 
Ref. || we found that we could select a subset of solutions of the equations of motion which 
is still a symplectic manifold. The method we used was to add by hand a set of second 
class constraints compatible with the equations of motion: this amounts to the selection of a 
symplectic submanifold of the symplectic manifold of solutions. We follow the same method 
here including the pseudo-classical spin variables 

As in the case of scalar particle, the Grassmann truncated semiclassical Lienard Wiechert 
solution A± s given in Eq.( |6.40| ) for the vector potential, with -k± s = E±s = — Jr^J-S for the 
canonical conjugate field momentum given by Eq. ( |6.45|) , provides us such a set of second 
class constraints by way of 

N 

Xi(r,a) = A±(t, a) -^QiA^ Si {a - ^(r), ^(r), &(r)) » 0, 

i=l 
N 

X 2 {t, a) = tu(t, a) - ^ Q l n ±Si (a - tfc(r), ^(r), &(r)) « 0. (6.48) 



i=l 



These constraints allow us to eliminate the canonical degrees of freedom of the radiation 
field and to get the symmetric semi-classical Lienard- Wiechert reduced phase space, in which 
there are only particle degrees of freedom. In analogy to what occurred in Ref. || for spinless 
particle the independent variables rfi, Ki and now & as well will no longer be canonical when 
one imposes these constraints by way of modified Dirac brackets. 

To determine the effects of these constraints in the construction of Dirac brackets we 
must compute the 6x6 matrix of brackets 



{Xi,Xi} {Xi,X2> 

{X2,Xl} {X2,X2} 



(6.49) 



It turns out that this matrix bracket is relatively simple, due to the Grassmann charges. 
Consider, for example the case of two particles. The particle or Lienard- Wiechert parts of 
the matrix bracket vanish since Q\ = = Q\ and cross terms vanish because they involve 
Poisson brackets of particle one variables with particle two variables. Thus the only part of 
the 6x6 matrix bracket that contributes is from the field variables. It has the form 

dd 

{Xi(r, e?i), X 2 (t, a 2 )} = (I - ^)5 3 (^ - a 2 ), (6.50) 

and since 

{%i,Xi} = = {%2,X 2 }, (6-51) 

only the 3x3 off diagonal portion contributes. 

In order to have a well defined Dirac bracket we need to use a modified form of the Poisson 
bracket in which the inverse of the matrix of constraint Poisson brackets is used. Calling 
this matrix C, we define C~ l so that CC^ 1 = (I — |f )5 3 (<xi — o^). But the transverse form 
of the delta function allows us to use the idempotent property of the projector to show that 
the inverse of C in this sense is just C itself. In that case for two functions /(/Cj, rfi),g(Ki, ffi) 
of the particle variables the Dirac bracket becomes 
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{f,gV = {f,g}- 

-Big[ J d 3 a{f, - E QiX ±Si {a - ff t (r), k^t), £>))} • 

i 

•{ _ 51 Qfi^Sji? - Vj(r), Kj(r),ii(r)),g} - 
j 

-{/, -EGi*W* - ^( r )' %( r )'^( r ))} • 
j 

■{- E - fh(r), iZi(r), ii(r)),g}] . (6.52) 

i 

This is the bracket for the new reduced symplectic manifold containing only particles. 
To find the new canonical basis for the particles with respect to these Dirac brackets, we 
define the following scalar function 

N-l N 

£ = E Yl QiQjfcij(€i,&,Ki,Kj',Vi - Vj), (6-53) 

i=l j=i+l 

in which there appear the following functions 
ICij = J d 3 a[A ±Si (a - ffi, &) • 7?^- (a - ffj, 

-^±Si(<? - 57j, Kj, £j) • TTlSi^ - ffi, K h &)] = Kijfe, £j) Ki, Kj\ ffi - ffj) = -JCji- 

(6.54) 

In Ref. we found that the old canonical variables ffi and need to be modified to 
ffi = ffi + T E QiQjdmfcij =Vi + Y, ^QiQj{ffi, fcij}, 

^i — Hi — QiQ jd^ JCij = + ~^QiQ j{^ii ^ij] 'j (6.55) 

so that they satisfy 

{ffk,ffi}* = 0, 
{kjs., = 0, 

{ffi, FZj}* = JSij. (6.56) 
In analogy to the above modifications we assume that 

i = d + \ E QiQilfi, £y}- (6-57) 

We now show that this definition, together with the previous ones, are sufficient to guarantee 
that the new variables are canonical. To this end we let both 

H"i l^i ~t~ 2 ^ / QiQj{[^ii 

Vi = Vi + - E QiQj{"i, JCij}, (6.58) 
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represent any of the new candidate canonical variables. Then we get 

{fii,Vj}* = 

- [ / £a{jl h - QkA ±S k(v - %(t), Kfc(r), |* fc (r))} 

•{- J2 Qrf±si(& - m( r ), £i( r )> 5( r ))> "A + 



d 3 a{fl h - QkK± S k(v ~ ff k {r), K k (r), &(r))} 

■{- H QzAlSf - ^( r )> 6( r ))> ^j} • 



(6.59) 



Since we have 



{f2i,Pj} = + ^{&,J2QjQk{Vj,fcjk}} - ^{Vj,J2QiQk{fii,K>ik}} + 

k^j - kj^i 



QiQkifii, fcik}, QjQdvji ^}}, (6.60) 

k+i l+j 

we see from the definition of /Qj that the last bracket would require that either % = j,i = I, 
or j = I. Due to the Grassman charges this forces the last bracket to be zero. For similar 
reasons we can replace the new canonical variables by the old ones in the last two lines of 
Eq.( |6.59 ) and with the sums truncating so that 

{fauj} = + ^{fii,Y,QjQk{pj^jk}} - ^{fy'lLQiQkiiUifcik}} 

k^j - k^i 

-QiQj J d 3 a[{[li,Axsi(v-Vi,Ki,£i)} ■ {n±sj(v -f}j,Kj,£j),Uj} 



n_LSi(v - Vi, Ki, &)} • {A_LSj{a - ffj, Kj, Pj} 



(6.61) 



Using the definition of the Kjk, we find that the second and third terms on the right 
hand side become 

k^j - k^i 

1 / r -> -» -» 



{fli, {z/j, 7i_L Si (a - ffi, K h &)}} • A_L k (a - ff k , K k , £ k ) 

{Vi, {fli, A ±Si (a - ffi, Ki, &)}} ■ 7r_L fc ((j - ff k , K k , Cfe) 



-{Vi, {fli, 7r_LSi(<? - ffi, Hi, &)}} ■ Aj_ fc ((J - ffk, K k , £ k ) 



+ 



+QjQi2 \ I da {^h^±i(^-Vk,^k,^k)}-{Pj,A±sj(^-Vj^j^j)}' 



-{f!i,A ± i(a -ffi,Ki,£i)} ■ {Vj,TT_LSj{v -Vj, 



-QjQ%2\ J d ^[{^j^^ji^ -Vj^j^j)} ■ {&, A ±Si (a - ffi, FZi, &)} - 
-{uj,A ±j (a -ffj,Kj,ij)} ■ 



(6.62) 
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Combining the terms proportional to 8ij by using the Jacobi identity and the fact that 
{z?i,/ij} has a zero bracket with both 7?j_si and A±si, implies that these terms vanish, while 
the remaining terms combine to exactly cancel the last two lines of Eq.( |6.6lD . Thus we 
have that 

= (6-63) 

and hence the new dynamical variables have the same Dirac brackets as the Poisson brackets 
of the old dynamical variables, that is, they are canonical. Thus, not only do we have same 
brackets as in Eq. ( |6.56P but also 

{/*,£•}* = {£,§} = 0. (6.64) 



As in Ref. || the rest frame condition 

N 



Up = V{int) =J2^i+ I d 3 a[7f ± x B)(t, a) = 
i-i J 

N 

= E ^ + E QiQj / d 3 a 7T± S i(a - ffi, Ki) x (d a x A_L Sj (a - ffj, %)) + 

i— 1 i<j 

+n± S j(a - ffj, x a x A±Si{v - Vi, « 0. (6.65) 

can be rewritten in these new canonical variables as 

~ N ~ 

ftp = V(i n t) = K+ = E ^ = °- ( 6 - 66 ) 

t=l 

This is due to the fact that the steps of expanding the cross products, integrating by parts 
and using the transverse gauge condition make no additonal reference to the spin dependence 
of the fields. 

The previous results allow to get V^ n t) in terms of the final canonical variables. For the 
internal angular momentum in terms of these variables we get 

N i N - - r 

JU) = e rst Sf = £(*(r) x ^(r)^ - - Ete x &r + / d 3 a (a x [ff ± xB](r, a)J = 

i=l z i=l J 

^ ~ ~ - 7 N % 

= YM - ck) x - A)] r - 5 Ette - li) x (6 - 7,)] r + 

i=l Z i=l 



1 * * r - 

+ 2 E E QiQj I d 3 a (a x [vr ±5j x x A ±5i )] (6.67) 

i=l j^i 



In Ref. we showed that —a* x — rfj x $ exactly compensates the last expression for 
the field portion (a, x $ turned out to be zero due to Grassman truncation). The arguement 
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we used depended explicitly on the fact that n± Si , A LS - have expected transformation prop- 
erties. In the case here of spinning particles, we need in addition the terms — 7^ x £3 — £j x 7$. 
Using the forms for d-j , (3 i7 and 7, given in Eqs.( |6.55| ), ( |6.57|) , together with the expression 
for K-ij, and expanding the cross products in the integral, using the transverse nature of the 
field together with vanishing surface terms, we find that 



N 



N 



J(int) = X R iY 



i=l 



-n£(&x&) p + 



N N 



+ 9 E E QiQj [(.Vi X d Vi +KiX 8 Ki + & X % 
A i=l #i 

cf 3 0-(Aj_Sj • 7T_LSj - • 7f_L 5i 

N N 



1 



E E QiQi / ^ 



^4±5« x n±Sj - 
-ax a A k ± si)ir k ±Sj + A ±Sj x n ±St + ax (A k xsAtt^ 



(6.68) 



In order to see how the cancellation between the last two lines works, we first spell out 
the new vector and scalar dependence of the vector potential and its canonical momentum. 
Their general forms are [from Eq. ( |6.40[) and Eq.( |6.45D 1 



.4 



±Si 



Kifi{ K i ) K i ' Pi) 4" Pi9i{ K i > K « * Pi) ' ^7i(Pi)^ni( K j j ' Pi) ~\~ 



, KiKi ' ' Pi . 1 2 -> \ 1 Pi^i ' 'Pi / 2 -> * 1 

H JiK , «i • Pi) H mi(«i , «i ■ pi, 



Pi 



Pi 



Arcp 



^[piCiirf, Ki ■ pi) + 



+ £S ■ Pi E e n(Pi)d n i(tf, Ki ■ pi) + £»& • «i E a n (Pi)b ni (Ki , «j • £j) 



(6.69) 



where 



ft 



1^ - ^ 



I a - ffi\ := 9 CT pi = —il-pipi 

Pi 



Now we use 

fji x (9^ = -ct x d a + pi x <9 CT , 
and several facts (using generic functions h): 

Pi x d a h(pi) = 



0. 



(Pi x ~\~ K>i x d Ki )hi(K%, ^ ■ pi) = 0, 

(pi x d a + KiX d Ki )Ri ■ = KiX, 



(6.70) 



(6.71) 
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(pi x d a + Hi x d K .)pi 

(Pi X d a + KiX d Ki + (i X • 
(Pi xi + KjX^+^X • ■ piKi 

(pi x d a + KiX d Ki +(jX d^)f>iKi ■ • p;/3i 
(pi xi + /j,x4+(iX ' 



PiX, 

K i ' ' Pi K x 5 

Pi^i ' Ki ' PiPjX. 

(-)£& ■ piii x • 



(6.72) 



where the x indicates vector cross product. Thus the effect of the bracket (pj x d a + k, \ x 
d Ki + x c^J is to turn scalar product into cross product, which together with the first 
part of Eq. (|6.68| ) leads to the cancellation referred to above yielding 



N 



N 



J{int) = J2(Vi X «i) - ly YJ& X 

i=l A i=l 



(6.73) 



Thus the total angular momentum of fields plus particles reduces to an expression in- 
volving just the new canonical particle variables when the fields are eliminated in pairs by 
using the modified Dirac brackets. 
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VII. THE SEMI-CLASSICAL HAMILTON! AN WITH THE DARWIN AND 

SPIN-DEPENDENT POTENTIALS. 



In this Section we study the Hamiltonian M of Eg. (|5.1|) , which contains both scalar 
(inside and outside the square roots) and vector (inside the square roots) direct interparticle 
potentials. By using the semi-classical property Q 2 = 0, these potentials can be reexpressed 
in terms of a unique scalar potential outside the square roots. This semi-classical potential 
contains the Coulomb potential and, moreover, semi-classical relativistic Darwin and spin- 
dependent potentials. The determination of these potentials will be done firstly by using 
the old variables (not canonical with respect to the Dirac brackets ( |6.52| )) and then in the 
final canonical variables. 



A. The Semi-Classical Hamiltonian in the Old Variables. 



If we rewrite the Hamiltonian M ( p.l[ ) in the following form 



N 



M = £ [ A 2 + (*i(r) ~ Qi Mr, Vi(rW + 

QiQj 



i=l 



+ 



+ §4tt | ffi(r) - ffjir) | 



■ {Vi{r) -%(r)) 



2 (rrii + Jmf + k 2 ) x mf + k 2 



QiQj 



N i 

:= \ m 2 + kAt) 2 + — . 



+ V DS (r 



(7.1) 



i=i 



with the scalar potential Vds(t) containing the generalized Darwin and spin- dependent 
interactions. This potential can be decomposed in three terms 

V ds {t) := U x (t) + V dso (t) + U(t), (7.2) 

where U(r), given in Eq. (|6.46 ), contains the contribution coming from the radiation field 
energy evaluated on the semi-classical Lienard-Wiechert solution. 
The potential U\(t) has the form 



tt ( ^ ^n i -2Q i H i (r)-l ± (r,ff i (r)) 
U\{t) = ========== . 

i=i 2Jm 2 + Kj(r) 2 



(7.3) 



containing the Grassman truncated vector potential, which was in the kinetic terms under 
the square roots, and the spin-magnetic field and spin-electric field interactions present in 
Eq.(PD 



Tii = iQS x & ■ B(t, ffi) - 2i 



rrii + J mf + k 2 



-2Q^(r)- B(r,fJ l (r)) + 



7r±(T,?7i(r)) x Ki(r) 



mi + dm 2 + k 2 (t) 



(7.4) 
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Finally, the potential Vdso{t) come from the interaction of the spin with the Coulomb 
electric field (after quantization it will contain important Darwin and spin orbit terms) 



Vdso(t) 



£ 



QiQj 



■ (Vi{i~) -%(t)) 



£ 



■£ 4tt I fk(r) - ffjir) | 3 ( m . + ^2 + + ^ 

QiQj ^i{r) ■ (ffi(T) - ffj(r)) x %(r) 



J 4vr I j£(t) - %(t) I 3 ( m . + Jm? + F%)Jm% + 



(7.5) 



The rest-frame conditions are given by the vanishing of the internal 3-momentum [the last 
line is Eq.( |6.66D 1 



N 



V(i n t) = H p (t) = k+(t) + / d 3 a[n ±s x B s ](t, a) = J2 R + ~ °> 



(7.6) 



i=l 



with the term coming from radiation field 3-momentum given in Eq.( |6.47| ). 

Using the closed forms for the vector potential and the electric and magnetic fields given 
in Eqs. fl6.40|) , ( |6.45 ) and ( |6.44 ), respectively, we can evaluate both the spin-independent 
and spin-dependent parts of the scalar potential U\ (r) of Eq. ( |7. 3| ) . 

While the last term in U\{r) gives 



i=l 
N 



-£ 



mf + /?j(r) 2 

QjQi 



& Ait J + Ki 



1 



\f)i - ffjUyJmj + Kj + ^m) + (kj -fjij) 2 ) 



the magnetic and electric terms in lZi(r) give 
+ x & ' B( T > Vi) _ 



m 2 + Kj(irij + Jrrij + k 2 ) 



■ [Kj + fjijKj ■ fjij] 
1 



(7.7) 



1 2\ mf + kAt) 



7v iQiQjii x ii ■ {[Kj i^ij ■ 8^+ j^^^'^^ +m ^ 



E 



xd,, 



1 



• j 1 



Vi ~ Vj\ Jrrij + (/% ■ fjij) 



nJm 2 + Ki 2 
?]}. 



(7.8) 



and 



A' 



-£ 



QSi ■ ■ n±(T,f}i) 



i=i \ mf + Ki 2 {mi + ymf + k. 
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N 



¥J \ m i + «t ( m i + \ m i + K 



6+ 



m 2 + Kj(rrij + J m ] + 



.Vii 



m 2 ^m 2 + K 



[m 2 + (/%•• ??ij) 2 ] 3/2 



i) - • ijij ■ dijKj ■ a. 



(7.9) 



respectively. 

In Appendix D there is a comparison between these terms for the potential U\(t) with 
the terms appearing in the potential U(r) of Eq. (|6.47 ). Introducing the differential operators 
kj, lj, fj and gj defined by 



k : — Ki + fj, fj - - 
li : = + gi, gi : : 



rrii + \ mf + Ki z )\ mf + « 



(mi + \/mf + Ki z )\/mi + k^ 



(7.10) 



and combining all terms, we find that the potential Vds(t) can be rewritten in the form 

Vds(t) = V DSO {r) + £/i(r) + C/(r) = K iD5 + £W (7.11) 



where 
Vlds 



E 



QiQj [" " "I - ^« " Vij^j ' Vij 1 



' 47T 



— z- 



K i * SIS? ^ij 



— « 



x 4"? ■ (kj x <9jj 



4y / 'm 2 + Ki 2 ^m 2 + Kj 2 Vij( m i + \/ m2 + ^i)\J m i + ^ 2 



1 



2ymf + Ki 2 Vij\J™j + (Zj-Vij) 2 ^m 2 + /^(m^ + ^m 2 + k?) 



Vij 



m 2 Jm 2 + K 2 



y/mj + i? j ( y /m? j + K} + m 3 ) Vij K + (*i • ^') 2 ] 3/2 
- - - - - 1 1 



-1) 



+ 



+ 



+ 



fj ■ ft + fi) 



^Jmf+k 7 / rjij ^mf + (Ki-fjij) 2 ym 2 + /% 2 ^m 2 + (% • r^) 2 

' fj "I - ' f i ~\~ f i ' f j 1 



2\/m 2 + K^ x jm 2 + Kj 2 %' 



+ 



+ \ Ki • fj + ■ f j + fj • fj (fj^j ' %j + fjAKj ' ?7ij ~l~ f/fi ' %7') ' %j 

( 



2 ym 2 + Ki 2 ^m 2 + Kj 2 r}ij ^m 2 + Kj 2 ( ^/m 2 + /^ 2 + ym 2 + (^ • r^-) 2 ) ' 
— 2((f6j ■ «i)(fj ■ /Cj - (fj • %-)(/«j • %)) - fj • [(/Cj • %•)(«, - %(/«j • %))](«« ■ %) + 
+ (fi • «i)(«j • Ki - (Kj ■ r)ij)(Ki ■ fjij)) - Kj ■ [(fi ■ flij)(Ki - f]ij(Ki ■ r)ij)))(Ki ■ fjij) + 
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+ (f; • ■ Ki - (fj ■ f}ij){Ki ■ fjij)) - fj ■ [(fi ■ fjij)(Ki ~ fjij^i ■ ' Vij)) 
1 1 ^ + 

VijyJ™?j + Z?\f™% + (^i • 2 (y/m^ + K. 2 + yjm% + • r)y) 2 ) 2 

+ • - Kj • (gj • (%%))] 



(7.12) 



The term Vlds above includes all of the terms in Appendix D that are in closed form 
and are of order 1/c 2 and higher while the term Uhds below includes all of the terms that 
involve the double infinite series and are of order 1/c 4 and higher. 

n.n. 00 00 r C. C. 
Uhds — 2^ 1^ 1^ / = ; ^ 2 ' / 2 ; ^ 2 

i^j 57F m=0n=0 ^mf + Ki 2 . yJlTllj + Kf 
K- K- 2n+2rn+l 
< ^) 2T " +1 ( 1 J = -di>) 2n+1 — 



m 2 + ^ 2 "Jmj + Kj* JJ (2n + 2m + 2)! 



kl . $..)( l 3 . d- )( = • d- ) 2m+1 

-* 2n+2m+3 
Kj q yn+1 'lij 



yJrrf + Kf J (2n + 2m + 4)! 

— * — * 

-( 11 • 4)( kj • 4)( & • 4) 2m+i 

^m 2 + Ki 2 ^jm^ + Kj 2 yJmi + Ki 2 

2n+2m+3 

/ 2j Q \2n+l "J I 

ymJ + S,- 2 ' ' (2n + 2m + 4)! 

'm 2 + /tj 2 \J m j + K j y m t + K i 



2n+2m+3 



K j q \2n+l 



^/m 2 + k 3 2 1 (2n + 2m + 4)! 



1 I — * — * 

K i K j / ^ \2m+2( K j S \2n+2 T hj 



2n+2m+3 



sjmj + Ri 2 ^m 2 + k 3 2 ^m 2 + ^ 2 W ^m 2 + * 3 (2n + 2m + 4)! 



ki * w k. 



, - 4)( " J = • <%)( = = • <%) 2m+2 

V m 2 + k, 2 ^ m 2 + Kj 2 ^jml + Ki 2 



2ra+2m+5 



"j 3 

As discussed in Appendix D this complex formula for the lower and higher order Darwin 
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and spin-dependent terms has been checked to reduce to the corresponding spinless results 
given in Eq.(6.14) in Ref. || when fj = gj = 0. 



B. The Semi-Classical Hamiltonian in the Final Canonical Variables. 



We now reexpress our final result in terms of the final canonical variables ( |6.56|) , (|6.64 ) 
forming a Darboux basis for the Dirac brackets. These are the physical particle variable in 
this reduced phase space. By using Eqs. (|6.55|) we get 



mf + 



\ 



mf + n,i + Ki • QiQjd^JCij 



Ki ^2jjLi QiQjdfjifcij 

mf + Ki H " 



(7.14) 



2\Jmi + Ki 



with 



QiQjK-ij^KiKjTji Vji^>ii^>j) QiQjfc-ij{Ki,Kj,Vi VjidiC. 



.))■ 



Ki 



dfjifCij = J d 3 a (pti ■ dfj.A± S i(Ki, ffi - fff, §)) • vr±5i(^, ?Zj, ffi - fff, £•) 



— A±Si( K i> K j> Vi ~ Vj'i £i) £j) ' ' ®iji'K±Si{Kii K ji Vi ~ Vji £j) 



(7.15) 



But as we have seen in Eqs. (|6.65| ) and (|6.66| ), we can replace J2i<j d Vi QiQjlCij by 
— J d 3 a[7Tj_s x Bs]{t, a). Using the integral in Eq. ( |6.47D , written in terms of canonical vari- 
ables, we obtain that the old kinetic term can be rewritten as the new final kinetic term 
plus an extra scalar potential U HDS (r/ij = \ffij\ = \ fji — %|) 



N 
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(2n + 2m + 4)! 



k -> 1 

' . d- )( 3 — • d - s 

2 ^2 ^ / ^2 ■ ?/ 

TO,; + Ki J m 2 _|_ ^ 



(( , "* " 4') 2m+2 ( 1 ' ' 4) 2n + 



to 2 + £ 2 



I / ^ a* \2m+l/ ^3 a* \2n+l\ ~2n+2m+3 , 



1 /m* I ^ yA / 9 ~2 °^V l / 2 , - 2 jJ 1 / m 2 7 -.2 °^ + 
V m i + K i Jm + k . V ™* + K t y m j + K j 



+ (-r==f • *) 2m (^=^ • ^) 2 " +2 >: 



~2n+2m+3 



V m i + K i ymj + k 



' L2 I r?2 ^ 1 / I ~2 ^ + 

to 2 + Ki V m * + K * y m] + K j 

kj ~ lj -*(^^-40 2m (^^-*) 2n+1 ])^ +2m+1 ]- (7-16) 



m 2 + /? 2 y m l + ^i . m 2 , g 2 



J 3 ' y 3 3 



There is some cancellations between the terms in U' hds (t) and those in Vhx>s(t) when 
expressed in the final variables. 

The final result for the semi-classical Hamiltonian is 



m = vi mt) = E ^mj + ^ + E — Ej f QiQAi&j + E + Vos(r) = 

i=i i=i 2ym 2 + Ki i^j Vi Vj I 



(int) 



= £ ^to 2 + * + E : , ~ 3 ~, + Vbs(r), [Vbs = Vbs + t^ DS ] (7.17) 
i=i 4vr | ?7i - | 

where 

Vds = Vlds + ^i/D5 + U' HDS := Vlds + Vhds, (7-18) 

with 



w=-E^EE[ 



^i) 2m+2 ( ) -^) 2 " + 
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Kj \2m( K j ft \2n+2\ ~2n+2rn+l 



yjrnf+rf* \jm 2 + K 



i 3 ---dlj) 



3 1 3 

— * — * 

/ Ki ft \2m+l( K 3 ft" \2n+l\ ~2n+2m+3 , 

^ j— — ^2 ' Ui i> \ I ' °' i hj 



m i + K i \\m) + K j 



j ^2 v/V / ~ * J / y V y v/ v / 



^2 

■3 



( ^ o*\2m+l/ \2ra+l\ ~2n+2m+3 , 

v j— =2 ' Ui i> y 1 Z2 ' 3 ' J 'v " + " 



+ 2( , 11 • 4)( ; = • 9, 



\/m 2 + Fz. 



m 2 + k^- 



_^ 5 n2to+1/ K j a" \2n+l~2rt+2m+3 

^2 ' Ui 3' \ / ' V' 'is 



( 2 ^ + 2 - + 2 ) !V v / STf 





~2 


\l m h 





k,; lj . Ki . n_ 1 1 , 



[ p = • ^(-7==* • ^-) 2m+1 ( i 3 = • + 

V "I? + Ki V m? + K^ 





~2 




~ K 3 



kj Kj 



-2 



V™F + «i \/m? + k- 



j 1 j 

kj ~ jj a* / a* \2m( S 0* \2n+lA ~2n+2m+l _ 

m 2 + k- ymf + Ki ^ m 2 + k- 

k- k- £■ 2n+2m+3 

xfm^t ' r^I^ xfm^t r7~^ (2n + 2m + 4)! + 



3 3 1 1 y 3 ' 3 
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£ 2n+2m+5 



Svf m+ \ , ' . ' ; 5Ti ■ (7-19) 



iriTf .L^ 5 " (2« + 2 m + 6)! 



The term Vlds is defined as in Eg. ( |7.12[ ) but with all the old non-canonical variables 
replaced with the new canonical ones. Due to Grassmann truncation the expressions are 
equivalent. As we shall show in the next Section, it reduces in lowest order to the pseudo- 
classical version of the Breit results (of order 1/c 2 ). The term Vhds involves the double 
infinite series and is of order 1/c 4 and higher. In general there is no closed form for it (see 
Section VI of Ref. [||] for related consideration for spinless charged particles), but in the case 
of the two body problem the double infinite series can be summed to closed form in the rest 
system. 
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VIII. THE SEMI-CLASSICAL 2-BODY PROBLEM: HYDROGEN ATOM, 

MUONIUM AND POSITRONIUM. 

A. The Semi-Classical 2-Body Problem. 

In this Section we specialize to the two body problem, which includes the pseudo- 
classical analogues of the hydrogen atom, muonium and positronium. 



1) The rest frame form of the energy in the unequal mass case (muonium) is (77 := 1 7/1—772 1 
and k, :— Ki — — re 2 ) 



M = \Jm\ + Jt + \jm 2 + ~/+ ®0 + V ds (t), (8.1) 

where Vds = Vlds + Vhds- From Eq.( |7.12|) we obtain for arbitrary masses 
Q1Q2 r re + (re • 77) 1 



Vi 



LDS 



47T 



2^m\ + &\jm\ + re 2 ^ 

« - 66 • ?7 . re • 66 ■ if 



— 7 , ; — I 



fj 3 (mi + V 77i 2 + re ) V m 2 + re 77 s (777,2 + ym% + il) \] m% + re 

.6 x 6 ■ (k 2 x <%) 1 .6 x 6 ■ (ki x £ ) : 

-1 ; — 1 ; + z — 



2ym 2 + re 2 rjy m 2 , + (k. ■ fj) 2 2ym2 + re 2 rjy m 2 + (k ■ r)) 2 

re • II 



m 2 + re 2 (mi + Y m 2 + re ) 



»c • 66 • 6fy j ^7 / raj y_ + k _ + 

1 1 y/mZ + fffa + y/nZ + 'tf) ^ K + (^) 2 ] 3/2 
-Nfi ■ 66 -d v K-d v - ' 



+7 



^ ^7772 + (« • 7}) 2 
« • 6 



m 2 , + re 2 (m 2 + \Zm| + re ) 



»g • 66 • 6fy j jj / ^1 V ^1 + « 2 

^m 2 + re* 2 (mi + V™i + ^ K + " OT /2 

+i|* 2 • 66 • d v K ■ d v \ — 1 = ) + 
^ y m§ + (re • z}) 2 



frk 2 1 , f 2 ki I 



\jm\ + re 2 77-y/ mf + (re • r}) 2 \Jm 2 + n 2 r\\jm\ + (re • t}) 2 
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At • f 2 — /C • f i + f i • f 2 1 



2y / m 2 + At 2 y / m 2 + At 27 ? 

+ (/t • f 2 - K ■ ft + ft • f 2 - (f 2 K • fj - ■ fj + f 2 fl • fj) ■ ?)) 

( 



2y m 2 + At 2 y m 2 . + At 2 ?? fjym 2 * + K 2 (ym 2 + k 2 + y m 2 + (/c • r}) 2 ) 
-2(/t 2 (f 2 • K - (f 2 • 7))(/Ci • 77)) - f 2 • [(/tj • 77) (k - t)(k • ?)))](«• r}) - 
-/t 2 (fi • k - (ft • 77) (k • 77)) - fi • [(« • 77) (« - 77(k • 77))] (k • 77) - 

-(fl • At)(/t • At - (K-fj)(K -fj))-K - [(fi • T])(k - l)(K ■ T)))](k ■ T)) + 
+ (f 2 • «)(« • K - (« •?))(«• ?))) - K • [(f 2 • ?))(« - • ?)))](« • 77) + 

+ (fl • «)(f 2 • K - (f 2 • 7))(K • 77)) - f 2 • [(fi • 77)(« - ?)(« • 7)))](« • 77)) 
1 1 ^ 

77^/ m| + K 2 \jm\ + (At • 17) 2 {\]m\ + At 2 + \jm\ + (k ■ fj) 2 ) 2 
~ 2 

+ [«-gi -K- (gi • (7777))]- 



) 



ml + At 



-[«•&-«• (g2 • (W>)] 



fjy m 2 + At ym 2 . + (k ■ fj) 2 y m\ + At + ym 2 . + (k ■ 77) s 
2 



\fm 2 + K 2 



fj\J m 2 -\- k \J m\ + (« • r/) 2 \J m\ + k + \j m\ + {R ■ fj) 2 



(8.2) 



In Appendix E we obtain the result 



87r(m 2 — ml) 



ki • k 2 - 

77 



(, 



r»2 + +k 



+ 



mf + k 



m,2 + k 



\ mj + K \ml + n \ml + (K-fj) 2 ^ m\ + (« • 77) 



mf + At 



+ 



TOj + At 

\ m 2 + At 



l)ki • h 



2 . - 2 _ _ _ _ 
^ - l)k 2 • h + lj • 1 2 



\ m 2 . + At 



~2 



1 / ym| + « 
9 V 



m| + (At ■ 77)" 



m 2 + At 



m 2 + (At • 77)' 



+((. 



m\ + At + ym 2 , + (/t • 17) 2 ym 2 + At + y m 2 + (At • fj) 2 

2 . ^ 2 ? ? - — 2 

mj + k _ lA kil 2 ■ ■7777At 

-2 ^ o + 



)+ 



\ mf + At 



mf + At k 2 • (li • fjfj)K li • (1 2 • fjfj)K 



\ m\ + At 



^2 



1) 



( 



(\Jm% + At 2 + y'ml + (At • r?) 2 ) 2 {\Jm\ + At 2 + \Jm\ + {R- fj) 2 ) 2 



+ 



) 



-) 
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•((. 



m 2 + k 
\ m\ + k 



l)[ki • k\ 2 ■ k - ki • (1 2 • {nfj + f)K))(K ■ 77)] + 



^ 2 - l)[k 2 • /cli • k - k 2 • (li • (£77 + • 77)] + 

\ ml + 



+ [r 2 • kIx • k - 1 2 • (fi • (nf) + t)k))(k • 57)]) 

1 



( 



m 2 + k 



\ m 2 2 + (k ■ r)) 2 (y/ m 2 + i + y^ 2 + (/? • r)) 2 ) 2 

1 



-(( 



m 2 + k 



-l)k 1 -(l 2 -r}r/)(/t-r/) i + 



+ 



— — ^ - l)ki • (1 2 • fjfj){K ■ f}) 2 + li • (1 2 • fjf})(R ■ fj) 2 ) 
\ m{ + k 



)(2 + 



m 2 2 + k + v/ml + (k • 17) 2 ) 2 ^ ™l + « 



m 2 + (« ■ fj) 2 



)(2 + 



(\/m 2 + K 2 + ym 2 + (/t-r/) 2 ) 2 \ mj + K 



m\ + (/c • f]) L 



7/ 



/-=, r ( k l - ll) • « - (kl — ll) • 77/C • 7) -, , - I K 2 - 1 2 ) • /V - ( K j - ! , ) • ///, // x 
— [K ■ K 2 K • Ki ~ 



(k 2 - 1 2 ) • K - (k 2 -l 2 ) ■ TjK ■ Tj- 



mf + K 



ml + k 2 



(. 



I \/ m 2 + K 2 yml + k 2 



\| m 2 + (k ■ f,) 2 y/ m 2 + / + ^ m 2 ~ ^ . ^ 2 



I ^2 A™ 2 I -»2 

y mf + K z ym2 + k z 



-ki • k, 



m\ + k + \Jm\ + (/c • 77) 
1 



\J 77Z 2 + (« -7/) 2 4- ^ 4- > -nr - ^ . nV^ 'I 



2\jm\ + n 2 \Jml + k 2 ?7 



m 2 + / - vm| + (k-77) 2 ) 2 (2 + 



m 2 + (« • fj) 2 



\ ml + k 



-Nm( + k 



V / m 2 + (/€-r}) 2 ) 2 (2 



+ 



m 2 + (k • 77) s 
\ m 2 + n 



))- 
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ki-k 2 ( 



m 2 , + k V m 2 , + k — y'm 2 .+ (R ■ fj) 2 



1712 + K + V m 2 + ' 5 7) 2 



m\ + R ^ym 2 + R - y m\ + (R ■ fj) 2 ^ ^ 

/ — > » ■ — — * — » — - — — — * — * — 2\ 

+ (^ki • £k 2 • k — ki • (k 2 • (£7) + f)R))(R ■ fj) - ki • (k 2 • 5717)/? J 



\jm 2 2 + (R-fj) 2 (\jm\ + + ym 2 , + (k • f)) 2 ) 2 
(^Jrn 2 + £ 2 ) 3 1 



mf + (k • 17) 2 (ym\ + R + Jmj + (R ■ fj) 2 ) 2 



— » — » • — o / 

-kik 2 • ■fjfj(R-fj) y 



)- 



m 2 , + R 2 (2y m 2 + R 2 + y m 2 + (/« • 77) 
(Vm 2 , + k + ym 2 , + {R-fj) 2 ) 2 



m\ + R 2 (2\Jm 2 + R 2 + y*m 2 + (R ■ fj) 2 



(y m 2 + R + ym 2 + (R ■ fj) 2 ) 2 



1 



mf + R^y'm^ + k 2 7? 



2) For the pseudo-classical version of positronium (mi = m 2 ) 



n/r n J 2 1 ^ 2 1 
M = 2y 77i + ft + 



47T17 



where now 



QiQ 2 [R + (R-fj) «-66 -?7 + «-66 -?7 



47T 



2(m 2 + n 2 ) 2 V r] 3^ m+ \J m 2 + /)\f m i + k 2 

1 



—2 



+7 



6 x 6 • (k 2 x <9 T; ) - 6 x 6 • (ki x ^) 



2\Jm 2 + k 2 

— # 



fj\Jm 2 + (k • 77) 2 



+ 



m 2 + K 2 (m + V 777 2 + k ) 
7K • 66 • 6^ 



fj m V 777 + K 



V / m 2 + /(m+ V / m 2 + /) ^ K + (^ " OT /2 



-1) + 



1 



+«6 • 66 • d v R ■ d n - 



V \Jm 2 + (R ■ fj) 2 



+7 



£•6 



(K 



\lm 2 + R 2 (m + \J m 2 + R ) 
iR ■ 66 • 6 4 



7/ 777 V 777 + K 



r77 2 + /(777 + V / 777 2 + /) ^ K + " OT 72 



1) + 



(8.3) 



(8.4) 
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'In 



V \Jm 2 + (k ■ fj) 2 



) + 



+ 



f r k 2 



1 



+ 



f 2 -ki 



2 + K 2 r]Jm 2 + (k ■ fj) 2 



+ k 2 rjJm 2 + (k ■ fj) 2 



K • f 2 - K • fi + fi • f 2 1 



2ym 2 + Kpym 2 + k 2 ^ 



+ [k • f 2 - k • fi + ft • f 2 - (f 2 K • 17 - fi/Cj • 77 + f 2 fi • 77) • 77] 

( 1 1 ) - 

2(m 2 + K, 2 )l] ^ m 2 _|_ ^.2^ m 2 _|_ £2 _|_ ^ m 2 _|_ . ^)2^ 

-2(/t 2 (f 2 • K - (f 2 • 7j) («j • fj)) - f 2 • • fj)(K ~ fj(K ■ fj))](K ■ fj) - 

-/(fl ■ K - (fx • ?))(« • 7/)) - fl • [(/C •?})(«- ?)(« • ?)))](« • 7/) - 

— (fl •«)(/€• K — (/C • 7)) (/C • 77)) — K ■ [(fi • fj)(K — fj(K ■ fj))](K • fj) + 

+ (f 2 • k)(k ■ K - (k ■ 7)) (« • 77)) - K ■ [(f 2 -fj)(K- fj{K ■ fj))](K ■ fj) + 

+ (?! • *)(f 2 • K - (f 2 • 77) • 77)) - f 2 • [(fx • 77) - fj(H • f}))](Z ■ fj)) 
1 1 

- + 



fjy 7n 2 + K 2 ym 2 + (k ■ fj) 2 (y m? + k 2 + \Jm? + [k ■ fj) 2 ) 2 

+ ([« • gi - K • (gi ' (W))] - [K • g 2 - K • (g 2 • 

2 . \jm 2 + tt 



fj\J m 2 + k Jm 2 + {k • fj) 2 \J m 2 + k + dm 2 + (/c • 77) 



8tt 



ki • k 2 - 



77 m 2 + (K-fj) 2 m 2 + - 



■) + 



+ (-fi-f 2 ; 



1 \l m 2 + k — \l m? + (k • 77) 



\J m? + k + \jm 2 + [R ■ fi) 2 \J m 2 + k \J m 2 + (k ■ fj) 2 



+ 



-)■ 



| li ■ (1 2 ■ ^ 1 

(\J m 2 + k + \jm 2 + {k ■ fj) 2 ) 2 \J m 2 + k \j m 2 + {k- fj) 2 

- (T 2 ■ «ll • K — h ■ (ll • (^T/ + 7}k) (« • 77) 

~2 

1 



m? + k 



\ m 2 + (K ■ f]) 2 ^ m 2 + ff + ^ m 2 + (£.^)2)2 



(- 



m 2 + k, J m 2 + (k • 7}) 2 



I ( L 
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lila • -fjfj(K-fjf 



-)(- 



m 2 + (k ■ fj) 



(\J m? + k + \j m 2 + (R ■ fj) 2 ) 2 \J m? + R \jm 2 + (R ■ fj) 2 
2 , 1 1 



m 2 + K m 2 + (R-fj) 2 m 2 + /" J r] 



- + 



+ ■ k 2 (li • « - (ki - fi) • fjR ■ fj) - k ■ ki(f 2 • k — (k 2 — I2) -fjK-fj) 
1 1 



( 



( 



\| m 2 + (k ■ fj) 2 ^ m 2 + £ + ^ m2 + (^-77)2 



1 

2fj + 



4(m 2 + k 2 ) 
1 



\ 



m 2 + (k • fj) 2 
^2 

m 2 + R 



( 



1 , 1 



m 2 + (R-fj) 2 m i + K 2 V 



A(\l m 2 + R — ym 2 + (is ■ fj) 2 )( 



2 + /t ym 2 + (R-fj) 



(ki-k 2 (- 



2(\Jm 2 . + R + J m 2 , + (R ■ fj) 2 ) 



(2\J m 2 + /5 — \/m 2 + (R • fj) 2 



m + k 



yAn 2 + (/? • fj) 2 



-(ki • /5k 2 • k - k x • (k 2 • (Rfj + fjR))(R ■ fj) - k x k 2 • -7777^ ) 

( . ... 1 ,. ) 



m? + k + \Jm? + (R ■ fj) 2 ) 2 



3y m 2 + R — 



-2 m + k 



(V m 2 + k ) ; 



^ m 2 + . 2(m 2 + (R ■ fjf 

ki • (k 2 • fjfj)(R ■ fj) 2i 



\jm 2 + R 2 



( 



(\/ m 2 + k + \jm 2 + (R ■ fj) 2 ) 2 



m 2 + (k • r)) 2 , 



m 2 _|_ £ Jm 2 + (R-fj) 2 2(m 2 + R ) (m 2 + R )fj 



(8.6) 



3) For the pseudo-classical version of hydrogen-like atoms (m 2 — > 00) the active part of 
the rest energy is just 



M=\lm 2 + ie + 9l9l + V 1 



4:711] 



LDS, 



1.7) 



since Vhds — ¥ in this limit. In this case Vlds simplifies dramatically to just 
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In the remaining part of this section we check the lowest order portion for comparison 
with the known Breit results. At lowest order (0(l/c 2 )) Eq.( [7.1 2| ) 



v lds (t) = J2 



QiQj 

Air 



K>i • Kj 4~ K i • fjij Kj • fjij 1 . K j ■ " Tjij 

I - 



i>£,i X £,i ' y^j X &ij 



Arriimj 



fi fi i 



+ 



2r] i jm i rrij 2m(nij rjij 

1 n 

rjijmimj - 
■ dij = i£i ■ dijii, Qi = so 



'2mimjT]ij 



+ 



+ ■ gj - K t ■ (cjj ■ {fjij fjij)) 



5.9) 



To lowest order f, 



V LDS (r) = y: 



QiQj r ^i ' ^j ~t~ ' Vij^j ' Vij 1 . ^i ' £,i£,i ' Tjij 



An 



Arriimj 



— i- 



Vij 



2rjf j m 2 



2r} i jm i m i 



+ & x 6 • (6 • d«f j x fy) 



" 2//,,///,/// 



+ 



+[-2mj • «f 2 6 ■ 4- + (& ■ ■ S)6 • 4- - 



(8.10) 

Including the derivative terms and canceling like terms we obtain our lowest order pseudo- 
classical expression, giving the semi-relativistic (order 1/c 2 ) potential energy terms for the 
pseudo-classical analogue of hydrogen-like systems (either mass — > 00), positronium-like sys- 
tems (mi = W2), or general unequal mass systems (muonium -like): 



V LDS (r) = J2 



QiQj 



Alt 



K j • Kj -\- Kj • fjij Kj • fjij 1 .Ki • ^i^i " Tjij 

7 ~ 



Arriimj 



Vij 



2ryfm 2 



13 1 



Xi X 6 ' X «j) + & X ^ ■ £ 2 X £ 2 1 + 47T^j) . 

2 1 1 1 j 1 / 1 j rjfj 3 



2rjfjmimj 



3& x ^ -£ 2 x r}y(^ • %)i 



2rjijmimj 



(8.11) 



B. Quantization of the Lowest Order Potential. 

Next we examine the quantum version of this interaction for comparison with the stan- 
dard results of the reduction of the Bethe-Salpeter equation. In the manifestly Lorentz 
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covariant formulation with the two signs of the energy the 5 Grassmann variables £^and £5 
are quantized to 7 5 7 M /v / 2and 7s/v / 2,so that the Dirac constraint — m£ 5 ps becomes 
the wave operator 75^7^ — m]. But now with the positive-energy spinning particles we have 
only 3 Grassmann variables £' r and no spinor wave equation (the spinors satisfy only the 
Klein-Gordon equation ). Therefore we have a different quantization (corresponding to the 
pseudo-classical spin S = (— z/2)£ x £ going into a/2) holding only on the positive-energy 
branch of the mass spectrum. The quantization is of Eq. (|6.64j ), {£f } = —iS hk , is 

r~^, (8.12) 

and is compatible with the Pauli matrix algebra of 

a h o k + o k a h = 25 hk . (8.13) 
With this conversion rules and with the algebra of the Pauli matrices we obtain 

{A- £i){B x g) = -hi xB)- % -A(B ■ Bi) + % -{A-B)a^ (8.14) 



giving 



Vlds{t) = £ 



4-k 

Vij X ' ®i _|_ Vij X K j ' &i 



l- 



AmimjT]ij Arjfjinf > ordered 



Arj^mj 2r)f j m i mj 
Bi-Bj ( \ t Airdiffij) ^ 3(Bi ■ Bj - B { ■ fjijBj ■ 



2mirrij rjfj 3 AmimjTjfj 



[8.15) 



For two particles in the CM system the above reduces to 
V LDS (r) 



Q1Q2 IfK 2 + (« • V) 2 1 . K-f) .K-f) 



4tt 



2m\m2ri rj 4r) 3 mf 4i] 3 my ordered 



L-Bi L- B 2 L ■ (Bi + B 2 ) 
4?7 3 mf 4r) 3 m2 277 3 m 1 m 2 

, 2vr ^ 1 B x ■ B 2 3 ^ ■ f)B 2 ■ V 1 1 , fi1R v 
+ "T^ Fi • °"2 + 7 5 7 5 • 8.16 

3 4 i] 6 4 m 1 m 2 J 

These terms, derived from the order 1/c 2 relativistic corrections beyond Coulomb poten- 
tial, are valid for hydrogen-like systems (either mass — > 00), equal mass (positronium-like 
systems) or general unequal mass (muonium-like systems). This expression is the same as 
derived from the Bethe-Salpeter equation and Breit equation (see Refs. |JT4]- |il||18| ) and by 



quantization of Wheeler- Feynman dynamics for particles with pseudo-classical spin flTi 
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IX. CONCLUSIONS. 



In this paper we have shown that the effective lowest order Hamiltonian used for the 
theory of relativistic bound states with spin 1/2 constituents of arbitrary mass (muonium-, 
hydrogen- and positronium-like systems) may be derived as a semi-classical result implied 
by the Lienard-Wiechert solution in presence of positive-energy spinning particles in the 
rest-frame instant form of dynamics. 

The result is non trivial in two respects: 

1) At least to us it was not evident that the semi-classical treatment of the delay with 
Grassmann-valued electric charges with Q\ = and QiQj ^ (regularizing the Coulomb 
energy and producing a unique semi-classical Lienard-Wiechert solution) could produce ex- 
actly a result derived from quantum field theory through the reduction of the Bethe-Salpeter 
equation. 

2) It confirms the validity of the concept of spinning particle, with its semi-classical 
description of spin with Grassmann variables, as a semi-classical simulation of fermions and 
the relevance of the Foldy-Wouthuysen transformation in determining the effective action- 
at-a-distance interparticle potential. 

Moreover, we obtain the full relativistic structure of the semi-classical approximation 
without any 1/c 2 expansions. This could help in evaluating relativistic recoil effects. 

We have already a semi-classical formulation for scalar quarks in terms of the Dirac 
observables in the non-Abelian radiation gauge in the rest-frame instant form of dynamics 
0. This has been obtained in the framework of weighted Sobolev spaces, where there is 
no Gribov ambiguity, in the case of a trivial SU(3)-principal bundle over each spacelike 
hypersurface. Having solved the Gauss laws, there is already a 1/r Yang-Mills transverse- 
potential- dependent potential in the final Hamiltonian. It acts between any pair of color 
charge densities, p a (r, o^) and p&(t, <t 2 ), [either quarks or localized charged Yang-Mills field 
(classical background of gluons)] and consists of three terms: 

i) a Coulomb potential between p a (j, <?i) and p&(r, CX2); 

ii) a potential in which p a (r, c?i) has a Coulomb interaction with an arbitrary color- 
carrying center located at 03, while p&(t, cr 2 ) interacts with the same center through a Wilson 
line along the flat geodesies 02 — 03 (plus the case 1^2): then one integrates over the 
location of the center; 

iii) an interaction with two centers <?3 and <?4 (over whose location one integrates): 
Pci(t, <7i) interacts with cr 3 through a Wilson line, the same happens between p&(t, <7<i) and 
04, and, finally, the two centers have a mutual Coulomb interaction. 

As shown in Ref. |7| in the case of mesons (quark- ant iquark system) the semiclassical 
approximation of Grassmann-valued color charges regularizes the potential and produces a 
semiclassical asymptotic freedom. However, due to the presence of the Wilson lines, it is not 
possible to check whether there is confinement. Neither baryons (3-quark system) nor the 
introduction of the spin along the lines of this paper have been studied till now. 

The challenge now is to see whether with the technology of this paper it is possible to treat 
the non-Abelian case of the quark model, by using a perturbative treatment based on the 
iterative Lienard-Wiechert-type solutions developed in Ref. [|33[ adapted to this non-Abelian 
radiation gauge. Namely, whether it is possible to get a non-Abelian Darwin potential and 
to check if it is a confining potential. 
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APPENDIX A: EXPLICIT FORMS OF THE F AND G FUNCTIONS OF 

SECTION II. 



To find the functions F — F 2 — -Fi, G — G 2 — G\ of Eqs. (|2.18| ) we must preliminarly 
determine the functions $ n , of Eqs. ( |2 . 1 2| ) , (|2.15| ). In order to calculate $ n and \l/ n from 
the definitions (|2.12| ),( p7L5| ), we observe that, since Eqs.( [2.6|) , (|2.7|) imply 



S = S free + 0(Q), 

S free = 2i(p ■ £)£ 5 6>(p), 

then, Eq.( p.l2|) implies that for every n we have 

2z£ $2n = {$2n-l&0, Sfree} => $2n = ~i(p • f)0(p)$2n-l, 
2*6^0^2n+l = {$2n£ , ^/ree} $2n+l = -{$2n, (p " 

Analogously for n > 2 Eq.( |2.15[) implies 

2i£ #2n = {*2n-l&£o, S 'free] ^2n = -l(p* • 0%) ^2n-l , 
2^56*2„+l = {^2n£0, %ee} ^2n+l = ~{^2n, (p ' 0%)}- 



(Al) 



(A2) 



(A3) 



From Eqs. (gTg) , (|2~T5| ) we get at order Q [{QA ^ ,S} = 2iQ$ 1 £ 5 £ , {p £o, S} 



*i - -'W • a? - 2 4r^ • Op • ^ + • a {-< h e^ J • (A4) 

Then, from {{QA £ , 5}, 5} = -4Q$ 2 &6,, {{p £ , 5}, S} = -AQ^ 2 ^ we obtain 
$ 2 = -^ 2 (p)(p • Of" <9A = f 2 , uv CC, 



m 2 = +i6 2 {p){p-i)i 



dA 

dx° 



(A5) 



where 



2, tin 



—i9 2 (p)p 



92,uv ■= +i0 (p)p 



dx° 



(A6) 
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1. Sum of Even Terms (Diagonal Terms). 

We prove the following lemma: 
Lemma 1 

Lemma 1 For n > 1 ; 

$2n = T {2n) (p){p-0£-dA {x), (A7) 

with 

Wp) = jmsr. m 

Proof: We first prove that, for n > 1 we have: 

$2n = f2n,uvCC- (A9) 



This equation is true for <3> 2 ; by assuming that it is true for <3>2n-2, then by using Eq. (|A2| 
we obtain: 



$2n-l = -0( P )p k f 2n-2,uvm V ,e} ~ {/ 2 n-2,™, RT^ = 

= ~id{p) [p U f2n-2,uvC ~ Cf2n-2,uvP V ) + K vk CC^, (A10) 

where 

h uvk = {f2n-2,uv,p k 0(p)}. (All) 

Using again Eq. (|A2|) and observing that £, u £, v £, k £, h = 0, we get 

$2n = -i9{p)(p-g)<$>2n-l = 

= -0 2 (P)(P- [p U f2n-2,u,C ~ Cf2n-2,uvP V ] ■ (A12) 

Thus Eq.flA|) is true with: 

f2n,uv = ~0 2 (P)P U [P k f2n-2,kv ~ f2n-2,vkp k ] , (A13) 

and thus by induction it is true for all n > 1. 

To complete the proof, we observe that the lemma is true for <3> 2 because: 



8 A 

f2, uv = T i2) (p)p^, (A14) 



with 

32 



T m {p) = -i6\p). (A15) 
Again by induction, since we have 
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f2n-2,uv = T(2n-2) (p)P 



dx v 



we can substitute this formula in Eg. (|A"9|) and, using Eq.( |A13| ), obtain 

$2„ = -0 2 (P)P U [P k f2n-2,kv ~ f2n-2,vkP k ] CC = 

= 9 2 (p)p 2 T (2n „ 2) (p)(p- £)£ • dA (x). 
Finally, by using Eq. ( |A15|) we get 

T(2n)VP) = (P)P T (2n _ 2) (p) = -% - . 



p- 



Q.E.D. 

In the same way we can to prove the following lemma: 
Lemma 2 



Lemma 2 For n > 1 : 



^2n = -T(2n){p)ip-m 



dA{x) 



dx° 



with: 



\p 2 e\p)f 

T(2n){P) = -« 2 • 



P 

The proof is the same of lemma 1, with $ 2n replaced by (n > 1) 

^2n = 92n,uvC 1 C , 

and f 2n , U v by 

,dA v (x) 



92n,uv = -T(2„)(p)p"- 



dx° 



Eq.( |A10| ) is replaced by 

*2n-l = -iO(p) \p U 92n-2,uvC ~ C<? 2 n-2,W>1 + KvkCC^, 

with 

Kvk = {g2n~2,uv,p k 0{p)}- 

Using the results of the two lemmas we find: 



^2n-^2n = -T {2n) ( P )(p-0i- 



dAjx) 
dx° 



+ dA (x) 



+T {2n) {p){p-i)i-E{x) 



Then, we have the explicit form of the F function given in Eqs. fl2.14| . pTT7| . |2TT 
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F — F2 — Fi — 

00 (2i) 2k 

= A (x)£ + (p- £)£ ■ E(x) £ -^T T (2k)(p) = 

• 00 o2fc 2kn2k( \ 

= A ( X )£ - & ■ E( X ) E(~) fc P l W 

P fe=1 //C! 

= A (x)£ - - 2 {v- Oi- E(x) [cos(2p0(p)) - 1] = 
= A (x)£ - i(p- 0(f • ^(x)) m ~/7" + _^ 
= A {x)£ o -i{p E{x)) 1 



m + + p 2 )^/™ 2 + 



2. Sum of Odd Terms (Skew-Diagonal Terms) 

We calculate first G x . In Eq.fl&lQl) for § 2n -i 

we can substitute the expression for fm,uv 

hn,uv = T {2n) {p)p U ^. 

Since we have 

h U vk — {f2n-2,uvi P k 9(p)} = h ukv) 

then, by symmetry considerations, we get 

K vk eee = 0. 

We have 

.( P 6(p)) 2n - 2 



—1- 



$2n-l = ~iO(p) 

{pO{p)f n - 1 ) 



p 2 



(p-dAo^p-O-p'i^dA, 



p 3 



p-dA ){p-£)-p 2 {£-dA )). 



Then we obtain 



00 „-2fc+l (n^atmWin-l 



G 1 = 2 i <p 1 + ((p.dA )(p.£)-p 2 (i.dA ))Y,( 5 ' i2l>lM ' l)y 



fc= -f '2k + l\ p* 
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' ' uy ^2A; + l! p 3 

, <!•/-> ^d6[jp)dAQ .2p6{p) -> -> 
+ 2 ^ • °^d^ + Z — {P -° {p - 9Ao) = 

= + i { i^± sin[2p0(p)] - ^ (gj^jjo) sin[2 ^ (p)] + 
L9 ._ xdOWdAo , 2pg(p) * 

Now we calculate G 2 . We have: 

and we can substitute the expression for g2 n ,uv into it 

d A v 

92n,uv = -T {2n) (p) P U — . 

With 

Kvk = {92n-2,uv,p k 0(p)}, 

we have in this case 

{dA v 1 
o^,p k e(p)h u ce = 

r) f)A v 



^3 



We get 



# 2n _! = +i0(p) 



_,(?*(P)) 

t 



2n-2 



p2 



i {pB{p)f^ xdFu h k 
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Then we have: 
G 2 = 2i*! + 

+ 



2A; + 1! 



. - 2 *)f • M - ^or - § + ^ - 6 ( + ^) + 

9a^j^ 1 p 2A; + 1! + 
dx° ^ p 2k + l\ 

= -e-^-H2^(p)] + 

■ 2pg(p) dA * i p,dF hk h k \ 

" Hr 1 ""'* ' « p ' a? + -S 2 !-^ ' ( +2ft V J ' (A38) 

The complete expression for the function G is 

G = G 2 - G 1 = -(!*• £)- sin[2p0(p)] + 

P 

+ ((p • ^)(p • a - \{p ■ a^^ fe ) ^ sm[2p^p)] + 

- ^U-Ap- 0(P- E) + ^|^. f) . (A39) 
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APPENDIX B: SPACELIKE HYPERSURFACES. 



1. Review of Their Properties. 

Let us first review some preliminary results from Refs. 0,^Oj needed in the description 
of physical systems on spacelike hypersurfaces. 

Let {E T } be a one-parameter family of spacelike hypersurfaces foliating Minkowski space- 
time M 4 and giving a 3+1 decomposition of it. At fixed r, let z^ir, a) be the coordinates 
of the points on S r in M A , {a} a system of coordinates on S r . If a A = (a T = r;a = {cr r }) 
[the notation A = (r, f) with f = 1, 2, 3 will be used; note that A = r and A — f — 1,2,3 
are Lorentz-scalar indices] and d A = d/da A , one can define the vierbeins 

4(r, a) = d A z»(r, a), d B z\ - d A z% = 0, (Bl) 

so that the metric on E r is 

9A6( T i °) = z a( t i °)V^z U b( t i 9tt(t, a) > 0, 



g(r, a) = -det \ \ g AB { r ^ a)\\ = (det \ | zj(r, B 



\2 



7(r, a) = -det\\g f s(r, a) ||. (B2) 



If 7 rs (r, a) is the inverse of the 3-metric g?§(T,a) [7 r *"(r, a)gi L s(r, a) = 51], the inverse 

J] 

7(r, 5) 



# AB (r, a) of ^(r, ^) [# AC ( r , °)9cb( T i °) = $i] is S iven b Y 



g Tf (r,*) = -[?-g T * 1 * f ](T,*), 
9 

g"(T, a) = Y s (r, a) + {lg Til g T ^ f ^]( T , a), (B3) 



so that 1 = g rC {j, a)gQ T {r, a) is equivalent to 
9(t, o) 



7(r, a) 



g TT {r, a) - ^ ts (t, a)g T f(r, a)g TS (r, a). (B4) 



We have 



and 



z?(t, 9) = (J ^> + ^7 fs "^)(r, <?), (B5) 
V 7 



^ = 4(r,a)^(r, ( x)4(r, ( T) = 

= (IT + z^)(^l, (B6) 



where 
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l 2 (r,a) = l, l,(r,a)z^r,a) = 0, (B7) 

is the unit (future pointing) normal to S r at 2 m (t, a). 

For the volume element in Minkowski spacetime we have 



<?z = z?{t, a)drd s ^ = dr[z^{r, a%(r, a)]J^{r, a)d s a 



g{r, a)drd 3 a. (B8) 



Let us remark that according to the geometrical approach of Ref. |34 



Eq.(p|) in the form zff(r,a) = N(t,B)V(t,B) + N f (r, o)z$(t, a), where N = = 
\J g TT — Y s 9Tf9rs and A^ r = g TS Y r are the standard lapse and shift functions, so that 
g TT = N 2 + g fS N f N 3 ,g Tf = g?sN\g^ = N~ 2 ,g^ = -N f /N 2 ,g fS = 7 « + 
£r = l^JL + Zs ,Y^, d±z = N^fdrd 3 a. 

The rest frame form of a timelike fourvector p^ is P M = rjy/p 1 ^; 0) = rj^rj^/p 2 , P 2 = p 2 , 

o 

where r\ = signp°. The standard Wigner boost transforming P M into p^ is 
JA,(P,P)=<£(«(P)) = 

P P - P + P 2 

= € + 2^(p)^(p) - ^ + ^))(«.(P) + 



1 + u°(p) 



v — e£(w(p)) = w^Qo) = p^/vyP 2 , 



v = r $(u(j>)) = {-uM-X- U ^ U &H (B9) 

o o 

The inverse of L M j,(p, p) is L fJ, u (P,p), the standard boost to the rest frame, defined by 

L^(p,p) = L^(p,p) = L» v (p,p)y^ f . (BIO) 
Therefore, we can define the following vierbeins Q 

eXp)) = J L^(p,P), 

<(«(P)) = ^(P.P)=^V^(«(P)). 

ej(«(p)) = V^o(u(p)) = u^ip), 

£(«(p)) = S r VrWp)) = (S rs u s (p); 5' - S r %h U "f^h , 

e}(v(p)) = u A (p), (Bll) 



27 The e£f(ii(p))'s are also called polarization vectors; the indices r, s will be used for A=l,2,3 and 
o for A=0. 
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which satisfy 



e>(p))£(«(p)) = 



4W P )y B {u(p)) = v i 



rT = e^u(p))v AB e u B Hp)) = <ApK(p) - £ e?(u(p)K(u(p)), 

r=l 

?7ab = e^(«(p))?7^e^(w(p)), 

Pa ^ (u(p)) = p< *^ (u(p)) = °- (B12) 

The Wigner rotation corresponding to the Lorentz transformation A is 
R^(A,p) = [L(P,p)A- 1 L(Ap,P)f I/ -- fi ° ^ 



R l j(Kp) 



J ""'" K " J /v(A ')"o • 



R l J (A,p) = (A- 1 Y 



f „a-1v {{A-')° -l)p p {A-y y 



p o + r] ^ A )°i Pp (A-i)P + ^ ] - (B13) 

The polarization vectors transform under the Poincare transformations (a, A) in the 
following way 

e?(u(Ap)) = (R-'yA^e'Mp)). (B14) 



2. General Form of the Constraints of Section IV on Arbitrary Spacelike 

Hypersurfaces. 

Let us try to undo the canonical reduction of Section III and to find the general form of the 
modified constraints ( |4.4| ) outside the radiation gauge and outside the Wigner hyperplanes. 
Since Eq. ( |3.35|) implies 



7T ± (T,a) = 7T [r,a) - — 



ar N 

i=i 



qt N N 

tT(t, a) + — QiP(? ~ Vi(r)) = vr r (r, a) + £ QiC r (r, 5), 

i=l i=l 



A = -& 



Ac(a) = S 3 (a), 



A 



5 3 (a) = c(a) 



47r|(j| 



d-c(a) = -5 3 (v) 



(B15) 
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and since Eq.( |3.39|) allows to write the Coulomb potential in the following form 

QiQj 



V. 



-QiQj / d 3 ac(a - rji(r)) ■ c(a - %(r)) 



-Q>Q.i I d'a^5\a - #(r)) • -^5\a - r^(r)) 



(B16) 



the general form of the modified constraint Ti! {j) ~ on the Wigner hyperplane without 
the restriction to the radiation gauge is 

n I ; 

n ' = e * " (E \j m l ~ iQi$(r)&{T)F r8 {T, ffi{r)) + (^(r) - Qj^T, ^(r))) 2 - 



QiQjKiir) ■ &(r) 



5 



E ^j'^w sn-y ^ ^(r) • ^ / d 3 <n c^i - ^(r)) • c(ai - fjj{r)) + 

¥j (mi + ymf + «j (r))y»7if + (r) 



(B17) 



Eq. ( |Bl(j| ) suggests the following modification of the ten constraints ( 3.16| ) describing the 
isolated system on arbitrary hyperplanes (namely not on the Wigner hyperplanes) 

H^t) = J dVft^r, a) = p S/1 - V{ 2 / ^{f^ °) + B 2 (t, a)) + 



A' 



+ V ™f - iQiCi{rMr)F rs {r, t*(t)) + +(^(r) - Q^ ± (r, ^(r))) 2 + 



»=i 



A' 



Oi(r)gi(r)^(r)^(r) 



^r 7 ^e 7 a (r)ef(r)6 T 



^ + N 2 ( r )\/ m I + ^ 2 ( T 
+ * L — / = = b T*hp{T)7r (r, ^(r)) + 



+ 



iY 



2 i= i ^mf + zt^r) 

TV 



■bu a {r)b^(r)Fm(T,f]i(T)) + 



<Oi(r)K«i(r)ef (r)6, a (r)ef (r)6 w (r) 7r'(r, *) + ? £f =1 Q fc (r)5 3 (a - ff k (r)) 



(mi + ^Jmf + K?(r))^/m? + £ 2 (r) 
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$ (r)b m (r) / d 3 a 



r d d 







r - N 

+ M r ){ / d 3 a[n x B] ¥ (r,a) + £ Mr) - Qi(r)Mr, %(?))]} w 0, 



W^(r) = 6jf(r) / d 6 aa r H v {r,a) - 6J|(r) / tfaWfaa) = 

= Sr(r)-(^(r)^-^(r)^){i|dW(7r 2 (r,a) + fi 2 (r,^ + 

1 AT 



i7rOiWOi(r)^(r)e?(r) 



m 2 + K 2 (r) A /m 2 + k 2 (t 



--b T7 Vs^f{r)^j{r)b Ta + 



N 



+ EW J + &( r ) - 0<(r)i4(T,^(r))]" + 



i=i 

N 



, i \/m 2 + /«f(r) 

^r ? r(r) gi( ; ) f (T ^ (r) 6,.,(r)?;, ;( r )/*',„.( r. ;/,( ~)H 



2 

i=i 



V m i + «i ( r ) 



tQiW^Mff (r)6, a (r)ef (r)^(r) 7T*(t, a) + f £f =1 Q*(r)5 3 (a - ff k (r)) 



gi(r)Q i (r)K tf (r)g(r)6, a (r) 
*#j (m i + x /m 2 + /? 2 ) x /m 2 + K 2 



ef (r)6 i/3 (r) / dW 



9 9 



5 



()a , A ^-m)\-^(a-ff J (r))} + 



+ mm{r) - bKr)bUr))\ I rfW[vf x B] s (r, a) + 



JV 



i=l 



0. 



(B18) 



However, there could be extra terms containing £|* b Tll [like the spin-spin and spin-electric 
field terms in Eq. (|3.16| )1, which vanish on the Wigner hyperplanes. To check whether they 
are needed, one should verify that the Poisson brackets of the ten modified constraints among 
themselves close on the constraints as it happens with Eq. (|3.16 ). We are not going to do 
this check here, since we are only interested in the rest-frame instant form of dynamics 

Again modulo terms in z tim (t, a) the modification of the original constraints of Eq.( j3.3| ) 
on arbitrary spacelike hypersurfaces is 



S3 



T-L^t, B) = p M (r, B) - l^r, B) 



=vr r (r, B)g rs (r, B)tt s (t, B) + 



+ 



2/y(r,<?) 
7 « (r> ?)7 «* (r> ? ) Fffi ( r> ^(r, a) + 



iV 



+ E*V-*0-))- 



i=i 



• V™ 2 - 7 r ~ s >, ^)(^r(r) - Q 1 {t)A,(t, a))(K iS (r) - Q^A^r, B)) + 
1 



N 



+ 



^(r)g j (r)e7(r)ef(r)/ 7 (r a)^(r)^(r)/ a (r, a) 



+ 



TV 



+ 7=,E^-^)) 



i=i 



Qi(r)ff(r)ef(r) 



l a {r, a)zsp{r, B)tt s (t, B) 



sjm\ - 7™(r, B)k^{t)k Cs (t) 

Qi{r)i?{r)^{r) 



-£E*V-*0-)) 
z i=i 



2 _ 7 rs( T) a)Ki?(T)Kis(T) 



■ Zua(r, B)z^(t, B)Y u (r, B)j sv (r, B)F h (t, B) + 
+ ^ 3 (a-^(r)) 

Qi(r)K if (r)^(r)zr a (r, B)^(r)z sp (r, B) [tt*(t, <t) + f £f = i Q*(^ 3 (<? - %(r)) 



dn r (-r)f?(T)2 fQ ,(r, B)^{t)z- sI3 (t, B) 



+ 



A 



+ J25\B- Ut)){kh ~ Qi(r)A,(r, B)) 
i=i 



(B19) 
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APPENDIX C: COMPUTATION OF THE FIELD ENERGY AND MOMENTUM 

INTEGRALS OF SECTION VI. 



Here we carry out the details in the computation of the field energy and momentum for 
the case N = 2. The general N results obtained in the text are an immediate generalization. 
From Eq.(fT33D and Eq. fTI^) we find that 



Q1Q2 



in' 



EE 



(2n)!(2m) 



■Vi-V 2 x 



Q1Q2 



m=0 n=0 

00 00 j 



EE 



8tt 2 

7 

(?h-4) 2m+1 



m =on=o (2n + 2)!(2m)! 



QiQ 2 



00 00 

EE 



1 



(V!-4)(U2-4) 

&-&)^V-fcl an+1 



^ 2 m =o„= (2m + 2)!(2n)! 



(V 2 -4)(Ui-^) 



(77 2 -4) 2n+1 )|a i -7f 2 | 



2n-l 



Q1Q2 

m=0 n=0 
2n+l> 



(2m + 2)!(2n + 2)! 



(^•9 CT ) 2m+1 )|a-r/ 1 | 2m+1 
(U 2 -4)(Ui-4) 



+ 



|2n+l 



9 CT (r/ 1 -a (7 ) 2m+1 )|a-r/ 1 | 



|2m+l 
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QiQ 2 

16tt 2 



1 



oo oo 

J? £(2n)!(2m + 2) 

I) 
2) 



V, x 



d a (\J 1 -d a )tf 1 -d a ) 2m+1 \a-f] 1 \ 2m+1 



+ (1 <- 



|2n-l 



+ 



(C3) 



Our aim here is to compute ~ / d 3 a(E]_ s + -B|)(r, a) and / <i 3 <7(-E_i_s x Bs)(t, a) . Essen- 
tially these integrals were computed in detail in Appendix B of || for scalar particles which 
correspond here to the result obtained by replacing V*, Uj — > 77^ . Following steps in that 
appendix similar to ones that lead to Eqs.(6.4) and (6.5) in |J gives the result Eq. ( p. 46] ), 
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APPENDIX D: COMPARISON OF U x {t) WITH U(t). 



Let us compare the terms ( |7.7| ), ( |7.8| ) and ( |7.9| ) of the potential U\(t) with the terms 
of the potential U{r) of Eq.( |6.46|) and let us try to combine them together so that, after 
the addition of Vosoij) of Eq. ([7.5|), we can find the expression of the Darwin and spin- 
dependent terms of the potential Vds(j) of Eq. (|7.2| ). We rewrite U{r) of Eq.( |6.46|) in the 
following form 



1JV (I D 00 00 /?• K • 2n+2m+l 

m =o„=o L ' K Jrf + W (2n + 2m + 2)! 
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mf + ^ J/ Jml + Zf 3 (2n + 2m + 4)! 



' V){ 3 t3){ W^k 2 " 3) (2n + 2m + 4)! + 



K i o \2m+l( K j % ^2n+l Vij 

I 2 + 

- - 2n+2m-l 

+ v, • v,( . Ki ■ 4) 2w ( . Kj ■ 4) 2n 7^ - 

' /^T^ \/m? + « J - 2 ' (2n + 2m)! 



2n+2m+3 



/?■ K - „2n+2m+l 

(V< • 4)(V,- • &•)( ; : • <9„Q 2m ( ; J = • di,) 2n -^ , (Dl) 

" ' ' ' lm 2 + A 2 \/m? + «i 2 (2n + 2m + 2)!-' 



3 1 J 



where we have rewritten the operators Vj of Eq.( 5.13| ) and U, of Eq. (|6.42|) in the form [kj 
and lj are operators] 



\Jm 2 + Ki 2 \jm 2j r Ki 2 (rrii + ym|~+£^)(m 2 + Ki 2 ) \Jm 2 + Ki 2 

U< = S + ^-4^-££ . = f ' ( D2 ) 

sJmf + Ki 2 (mi + ^ m 2 + Ki 2 )(m 2 + Ki 2 ) ^m 2 + Ki 2 

We consider the last two lines in the expression for U(r) given in Eq. ( |6.46| ) 



1..N 



QiQj \^ \^ r kj kj 
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ki - kj 



m 2 + Ki 2 Jm 2 + K j 2 
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mf + Ki 2 



■a, 



\ 2m i 



K J 



71 



QiQ 

+ ( 



\/ m i + >% 2 
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•a. 
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•a 

ki 
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(2m)! 
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• d l3 ) 2m + 
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2ml nj 
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kj 
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4)(E ([( 



m=0 
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■ d l3 y m + 
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AW 
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kj 
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\ 2n+2 'lij 



(2n + 2m + 6)!' 



(D3) 



The first three lines on the right hand side contains two separate types of summations, 
each of which can be summed to close forms. The first type involves 



E 



K ; 



m=0 Jmf + Ki 2 



. Q \2m r Jij_ 



2m-l 



1 



m? + Rj 2 



( 2m ) ! %' ^/mf + (k\ ■ fjij)' 
In order to perform the second type we use 



(D4) 



dij dij 



2m+l 



mf + K^^mf + Ri 2 (2m + 2)! 

(2m + l)!!(2m - 1)!! 
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-(i[K 2 -(^-^,) 2 r + 



r)ij(2m+ 2)!( v /m 2 + k 
+2m[KiKi - (Ri ■ + (k 2 - (Ri ■ rjij) 2 )" 1 ^ ~ 

-(l - 5 m o)mjVij(K 2 - (Ri ■ Vij) 2 ) 171 ' 1 ^ 2 ~ ( 2m + • Vij) 2 ))- 



(D5) 



where I is the unit dyad. This complex form was obtained by examples for m = 0, 1, 2, 3, 4. 
As a simple check on this result consider (for m ^ 0) the contraction Eq.( |D5| ) with the 
dyad RiRj . This gives 



SS 



mf + Ki 



2m+l 



(2m + 2)! 



(2m + l)!!(2m - 1)!! 
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(D6) 



and it agrees with the equation described above Eq. ( |6.32| ). Using Eq. (|D5|) , the second sum 



(we separate out the m = term and let x 



is 



E 

m=l 
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has three parts, the first of which 
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while the second part is 

" (2m+ l)!!(2m - 1)!! 
m =i %-(2m + 2)!( v /m 2 + ^ 2 )2™+i 

2m[/«ft - ft • rjij)(KiVij + Vij^i)} <A 2 - ft • ^j) 2 ) 



(D7) 



2\m-l 



- ft ■ Vij)(^iVij + VijZi)] g (2m + l)!!(2m - 1 ) !! 2 mx m 



ft ' 



mf + m=l 



(2m + 2)! 



(D8) 



89 



Using 
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Eq. ( |D8| ) becomes 
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Finally our third part of the sum is 
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Using the results above we get 
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Combining these three parts of U (r) gives us 
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Note that the double dyadic dot products are defined such that the right most k's are 
contracted with the right most fj's. 

Now we combine this with the expression (7/7) for the corresponding portion of U\ 
(involving A±). First note that 

_ j2 Qj^i( T ) ' A±(t, ffi(r)) _ 
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Now combine Eq.( pi7| ) with the above term Eq.( pi5[ ) from U(t). In order to clarify the 
cancellations that will take place consider the case without spin for which the spin operator 
containing terms kj and lj reduce to the simple momentum factor Kj. In that case Eq.( pi5 ) 
becomes 
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while the corresponding portion of U\(t) becomes 



QjQ 



i+j 



71 



Kj ■ Kj 



yjm% + k 2 rjij ^m 2 + (Kj • fj^y 
1 



+ 



yfrrfi + k- 2 Vij yfrn% + (/* • %) 

[Ki ' K j K j ' V'ij^j ' Vij] 



rjij yJrn% + K? yj 'm 2 + (kj ■ 



-[Hi • Kj • fjijKj ' ^jj] 



^m 2 + k 2 yjm 2 + • 77ij) J 



\Jm 2 + K 2 + yJm^ + iKi-Vijy 



Adding the two gives 
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(D20) 



which is our relativistic extension of the standard Darwin interaction for spinless particles 

Now that we have checked our results without spin, we combine the spin-dependent terms 
of Eqs. (PI5Q , ( FT7D to get (k := k + f;l:= k + g) 
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To this we must add the remaining parts of U\(t) and U(t). To complete U\(t) we need 
the expression for IZi in Eq.( |7.4| ). The magnetic portion (|7.8| ) is 
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while the transverse electric part ( |7.9j ) is 
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(D24) 



Combining Eq.(|D2g) and Eq.(p24|) with Eq. (|D2ll) , with Eq.(^D for V DSO and with the 
first four lines of Eq.(|Dl|) produces Eq. ( |7.11| ). 
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APPENDIX E: SUMMATION OF TWO-BODY REST ENERGY TO CLOSED 

FORM 



In the special case of the two body system we can obtain a closed form if we use the 
rest frame condition K\ + k 2 = . The expression we get in this way may be used with the 
DiracjDrackets associated with R\ + k 2 ~ 0, so that the final reduced phase contains only 
fj — I ?7i_ — 77b I and k := K\ = —k 2 . Using the identity below Eq.( |6.32| ) 
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and the identity Eq. (|D5|) in the form 
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the higher order Darwin plus spin part Vhds — Uhds + U' HDS becomes 
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and similar identities applied in Appendix D and we obtain Eq. Q8.3p . 
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